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Abstract 

This work summarizes and puts in an overall perspective studies done within CMS concerning 
the discovery potential for squarks and gluinos, sleptons, charginos and neutralinos, SUSY 
dark matter, lightest Higgs, sparticle mass determination methods and the detector design 
optimisation in view of SUSY searches. It represents the status of our understanding of these 
subjects as of Summer 1997. 

As a benchmark model we used the minimal supergravity-inspired supersymmetric standard 
model (mSUGRA) with a stable LSP. Discovery of supersymmetry at the LHC should be rel- 
atively straightforward. It may occur through the observation of a large excesses of events in 
missing Et plus jets, or with one or more isolated leptons. An excess of trilepton events or of 
isolated dileptons with missing Et, exhibiting a characteristic signature in the invariant 
mass distribution could also be the first manifestation of SUSY production. Squarks and gluinos 
can be discovered for masses in excess of 2 TeV. Charginos and neutralinos can be discovered 
from an excess of events in dilepton or trilepton final states. Inclusive searches can give early 
indications from their copious production in squark and gluino cascade decays. Indirect evi- 
dence for sleptons can be obtained also from inclusive dilepton studies. Isolation requirements 
and a jet veto would allow detection of both, the direct chargino/neutralino production and of 
directly-produced sleptons. Squark and gluino production may also represent a copious source 
of Higgs bosons through cascade decays. The lightest SUSY Higgs h — > 56 may be reconstructed 
with a signal/background ratio of order 1 thanks to hard cuts on justified by escaping 

LSP's. The lightest supersymmetric particle of SUSY models with conserved i2-parity represents 
a very good candidate for the cosmological dark matter. The region of parameter space where 
this is true is well-covered by our searches, at least for tan/? = 2. 

If supersymmetry exists at electroweak scale it could hardly escape detection in CMS, and 
the study of supersymmetry will form a central part of our physics program. 
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1 Introduction 



There are strong arguments that the Standard Model (SM), despite its phenomenological 
successes [0, is only a low energy effective theory of spin-1/2 matter fermions interacting 
via spin-1 gauge bosons. A good candidate for "new physics" is the Supersymmetric 
(SUSY) extension of the SM which in its minimal version (MSSM) doubles the number 
of known particles, introducing scalar and fermion partners to ordinary fermions and 
bosons and relating their couplings. This scheme provides the necessary cancellation of 
quadratic divergences which appear in loop corrections to the masses provided the masses 
of super-partners are of the order of electroweak (EW) scale, i.e. ~0.1 to 1 TeV. 

One of the main motivations of the experiments at the Large Hadron Collider (LHC) 
is to search for SUSY particles ("sparticles"). Previous studies [0] have shown that at 
the LHC the observation of an excess in event rates of specific final states over the SM 
expectations would signal the production of SUSY particles. The present study is based 
on a more precise description of the Compact Muon Solenoid (CMS) detector and covers a 
wider spectrum of physics channels. The primary goal is to understand the experimental 
limiting factors for SUSY studies and contribute to the detector optimisation while its 
design is not entirely frozen. The emphasis is put on leptonic channels, which are as 
interesting in terms of discovery potential, if not more, as the classical multi-jets + 
^miss signature. The latter is more sensitive to instrumental backgrounds which can be 
fully understood only once the cracks, dead areas and volumes of diminished instrumental 
response due to detector services, mechanical supports are fully specified and evaluated. 
Another motivation for this study was the need to evaluate the discovery potential of an 
LHC detector in terms of accessible sparticle spectrum, sparticle mass reach, possibilities 
to determine the SUSY model scenario at work and possibilities to determine sparticle 
masses and model parameters in comparison and competition with a future e^e~ collider 
and dedicated SUSY (WIMP) dark matter search experiments 

This work summarizes and puts in an overall perspective specific studies done within 
CMS concerning the observability and discovery potential for squarks and gluinos 0, ^, 
sleptons 1^, charginos and neutralinos |^ H, SUSY dark matter [|1^], the lightest 



Higgs [Q, sparticle mass determination methods |10|, 0, detector design optimisation 



in view of SUSY searches [0, etc., and represents the status of our understanding of 
these subjects as of summer 1997. 

We discuss the specific SUSY model employed in section 2 and the experimental 
signatures investigated in section 3. Detector issues and simulations of detector response 
are described in section 4. Searches for SUSY in lepton{s) +jets + i?™** final states from 
strongly interacting sparticle production are discussed in section 5 followed by inclusive 
searches for the next-to-lightest neutralino in section 6. Exclusive 2 leptons + no jets+ 
j^miss exclusive 3 leptons + no jets+ [EJp'^^''^) channels are discussed in sections 7 and 
8. Possibilities to observe h ^ 66 in squark and gluino decays are discussed in section 9. 
The methods to measure sparticle masses and restrict model parameters are addressed in 
section 10. Results and conclusions are summarized in section 11. 

The potential of CMS to study the Higgs sector of SUSY has been extensively discussed 
elsewhere [|14| and we do not address this problem here. Obtained results are based on 
two-loop level calculations |T^. To summarize, in the search of the various MSSM Higgs 
bosons, most of the m^, tan/3 parameter space will be explored through at least one 
channel, provided an integrated luminosity of Lint = 10^ pb~^ is available. A difficulty 
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still exists in the region of 120 GeV < < 240 GeV and 2 < tan(3 < 10, which would 
require L^nt = 6 -^ 10 x 10^ pb~^ to be fully explored (though section 9 indicates that the 
Higgs sector may also be detectable in the decays of other SUSY particles). 



2 SUSY model employed 

Since the phenomenological implications of SUSY are model-dependent, the discovery po- 
tential of a detector has to be studied resorting to some particular model, preferably with 
a limited number of free parameters. This implies some loss of "generality" , but ensures 
tractable predictions. Once the possibilities and problems have been well understood 
within a defined scheme it is easier to evaluate its generality. 

We have chosen to make the investigation in the minimal Super Gravity constrained 
version of the MSSM (mSUGRA) pl as implemented in ISA JET 0. This mSUGRA 



scheme has a limited number of parameters, has well established and simple mass re- 
lations, and is implemented in event generators, besides having definite phenomenolog- 
ical/theoretical attractiveness and plausibility. In the mSUGRA model only five extra 
parameters, which are not present in the SM, need to be specified: the universal scalar 
mo and gaugino mi/2 masses, the SUSY breaking universal trilinear coupling Aq, the ratio 
of the vacuum expectation values of the Higgs fields tanP and the sign of the Higgsino 
mixing parameter sign{ii). The first three parameters are fixed at the gauge coupling 
unification scale and the sparticle spectrum at the EW scale are then obtained via renor- 
malization group equations |ll8 |. 



Throughout this study except for cases specially mentioned we largely limit ourselves 
to the set: 

tan/3 = 2, Aq = 0, < 0, (1) 

and also consider five representative parameter space points suggested by theorists, and 
listed in Table 1.1. 



Table 1.1: Representative mSUGRA points suggested by theorists. 



Point 


mo (GeV) 


mi/2 (GeV) 


Ao (GeV) 


tan/5 




1 


100 


300 


300 


2.1 


> 


2 


400 


400 





2 


> 


3 


400 


400 





10 


> 


4 


200 


100 





2 


< 


5 


800 


200 





10 


> 



In the following we discuss the mass relations among sparticles in mSUGRA, display 
isomass curves in terms of model parameters, discuss production mechanisms and decay 
modes accessible to experimental observations. 



2.1 Charginos and neutralinos 

The mixing of the fermionic partners of the electroweak gauge and Higgs bosons, the 
gauginos and the Higgsinos, gives rise to the physical mass eigenstates called the charginos 
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(Xi 2) ^^^d neutralinos (x^ 2 3 4). They arc among the hghtcst expected SUSY particles and 
therefore present particular interest. The two lightest neutralinos and the lightest chargino 
(Xi) X2' xf) have as their largest mixing component the gauginos, and their masses are 
determined by the common gaugino mass, mi/2- Within mSUGRA: 



0.45mi/2 (2) 
iM^o (3) 
(0.25 ^ 0.35)Mg (4) 

Figure 2.1 shows the isomass contours of {xi, xf) and of (xg, x1, xt) in the 
(mo, mi/2) plane. We remind the reader that in most of parameter space and in the more 
plausible scenarios the lightest neutralino, Xi) is the Lightest Supersimmetric Particle 
(LSP) and as such is particularly interesting as it is one of the most plausible cosmic dark 
matter candidates particle physics provides. 

The lightest chargino Xi has several leptonic decay modes, giving an isolated lepton 
and missing energy due to the undetectable neutrino and LSP (x? in mSUGRA): 



xf 



XU 



three-body decay 



•xt 



•xf 



xt 



It ^ 

-X? 
^X> 



X? 



two-body decays 



Figure 2.2 shows the branching ratios for the above listed xf decay channels as a 
function of (mo,mi/2). One can see that in different regions these different leptonic 
decay modes are complementary, amounting to a larger than 20% branching in the entire 
{mo, 1711/2) plane (for I — e, /i). Three-body decays are open in the region mi/2 ^ 200 
GeV and mi/2 ^ 0.5mo, whereas in the rest of parameter space the two-body decays are 
dominant. 

Leptonic decays of X2 give two isolated leptons and missing energy (xi): 



X'2 



•X'2 



X? 1- 



xt l""^ 

7± 



xli 



V 



three-body decays 



X^ 



X? 



two-body decay 



Figure 2.3 shows the X2 leptonic branching ratios as a function of (mo,mi/2)- Like 
the Xa case, the three-body decay branching ratios are sizable at mi/2 ^ 200 GeV and 
^^1/2 ^ 0.5mo. Beyond mi/2 — 200 GeV, the three-body decay of X2 is suppressed due to 
the opening up of the channels X2~^Xi h and X2^ Xi Z ("spoiler" modes). The two-body 
decay branchings are significant if mo ^ 0.5mi/2. 
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In hadronic collisions, the charginos and neutralinos can be produced directly via 
a Drell-Yan mechanism or, more abundantly, through the cascade decays of strongly 
interacting sparticles. There is also the possibihty of associated production xf/Xj + 
g/q. There are 21 different reactions (8 xf'x'jj 3 xfxf and 10 Xi'x^) for direct chargino- 
neutralino pair production among which X2 production followed by leptonic decays 
of xf and X2 is most easy from the detection point of view, yielding a 3 leptons + 
no jets+EJp'^'^^ event topology. 

Figure 2.4 shows the cross-section times branching ratio, ax B(3/ + invisible)^ where 
a is the XiX^ direct production cross-section and B(3/ -|- invisible) is the convolution of 
all the leptonic decays of xt and X2 listed above. Within a relatively small region of mi/2 
and for mo ^ 400 GeV, the crx B(3Z + invisible) drops by several orders of magnitude 
and vanishes at mi/2 — 200 GeV. For < 400 GeV the slope is less sharp thanks to the 
presence of two-body decay modes. 

Probabilities for production of Xi and X2 from gluinos and squarks are shown in 
Figs. 2.5a and 2.5b for tan/? = 2, < 0. One sees that these two sources of charginos and 
neutralinos are complementary in (mo,mi/2) parameter space. This abundant production 
of X2 from strongly interacting sparticles over the whole (^0,7711/2) plane has very useful 
experimental implications, which will be discussed in detail later in this paper. The 
convolution of the X2 indirect production cross-section with its leptonic decay branching 
ratio is shown in Fig. 2.5d. Clearly, over a large portion of the (^0,^1/2) plane, one 
expects rather large rates of same-flavor opposite-sign dileptons originating from xH- 

Similarly to the X21 lightest chargino xf is predominantly produced from strongly 
interacting sparticles, as illustrated in Figs. 2.6a,b. Moreover, the leptonic decay branch- 
ing of the Xi always exceeds 0.1 per lepton flavor (Fig. 2.6c). For low values of mo, in 
the region where decay to sleptons are kinematically allowed, the Xi decays to a lepton 
with a probability close to 1. 

Figure 2.7 shows the sum of xH ~^ Xi^~^^~ and X2 ^ ^ Xi^~^^~ branching ratios for 
larger values of tan/3 = 10 and 35 for both signs of fi. These branchings to e"'"e~ or fi'^n~ 
are smaller than at tan/3 = 2 shown in Fig. 2.5c. 

2.2 Sleptons 

The SUSY partners of ordinary leptons are scalars. Left and right-handed charged slep- 
tons are not mass-degenerate. The slepton masses are determined by mo and mi/2 0: 

mf^ =ml + 0.15m^/2 - sin'^6wMlcos2f3 (5) 
mf^ =ml + 0.52m?/2 - 1/2(1 - 2sin'^ew)Mlcos2f5 (6) 

m? =mo + 0.52mi/2 + l/2M|cos2/5 (7) 

Asymptotically, for mi/2 —>■ the slepton mass is determined by mo. Charged left 
sleptons are the heaviest. The mass dependence on mi/2 for left (L) sleptons is stronger 
than for right ones (R), hence the left-right slepton mass splitting increases with mi/2. 
Figure 2.8a shows slepton isomass contours in the (mo,mi/2) plane. 

Since xf, X2 isomass contours (see Fig. 2.1) are approximately straight lines of constant 
mi/2, as one moves up a slepton isomass curve, the X2 niass increases and eventually 
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becomes higher than that of the slepton. Thus we may identify two domains in the 
(mo,mi/2) parameter plane: one is where left sleptons are more massive than xf, X2 
(mo ^ 0.45 • mi/2, domain I), and the other, where they are lighter (mo ^ 0.45 • mi/2, 
domain II). These mass relations, shown Fig. 2.8b are responsible for distinct contributions 
to slepton production and decay mechanisms. In domain I sleptons can be produced 
only through a Drell-Yan mechanism (direct slepton production), via qq annihilation 
with neutral or charged gauge boson exchange in the s-channel. Pairs with the same 
"handedness" can be produced, namely, Il^l, IrIri ^^^i^^iIl- In this part of parameter 
space both direct and cascade decays of sleptons to LSP are possible: 

i) the left sleptons can decay to charginos and neutralinos via the following decays: 

^ ' ^ Al,2 
It ^^I + Xf 

ii) for right sleptons, only decays to neutralinos are possible; they dominantly decay 
directly to the LSP: 

iR^i- + x'i 

If decays to the second neutralino or lightest chargino are kinematically allowed, the 
decay to the LSP can proceed through several decay paths (cascade decays) with missing 
energy and/or isolated leptons in final state: 

x° - X? + i^i- 

x^-x? + ^° 
xf^xl + i^ + i' 
xt^x'i + w^ 

In domain II indirect slepton production, from directly or indirectly produced 
charginos/neutralinos is also possible: 

x'2 lij^ 
X2 ^i^i 

Xt - i>il^ 
xt - lU 

When allowed, indirect I production arising from decays of strongly produced g, q to X2 
or xf, is the predominant production mode. With increasing mi/2 along a slepton isomass 
curve, the fraction of direct decays into the LSP increases as the mass difference between 
sleptons and charginos-neutralinos decreases, until in domain II sleptons can only directly 
decay to the LSP. Sneutrinos in this domain thus decay totally invisibly, i.e. no component 
(z/, x5) from sneutrino decay can be directly detected. The slepton production and decay 
features described above give rise to several characteristic experimental signatures: 

• two leptons + i^™'^* + no jets event topology arises from direct charged slepton pair 
IlIl-i IrIr production, followed by direct decays to the LSP. The same signature is 
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found in production when one of the sneutrinos decays invisibly and the other 
has a cascade decay with two leptons in final state. This signature can also be 
found when sleptons are produced among decay products of directly or indirectly 
produced charginos-neutralinos. In the case of indirectly produced sleptons there 
are the following differences compared with direct production: i) single production 
is also possible, for example in production followed by decay of the ^ 

slepton; ii) a Dili pair can be produced from direct xfxi"; iii) left handed sleptons 
production is favored; 

• three leptons + i?™'*^ + no jets final states are expected only from vil^ production 
followed by cascade decays to LSP; 

• four leptons + Ej''""" + no jets comes from IlIl-, ^i^-, ^iIl production followed by 
cascade decays to LSP. 

The last two topologies can also be found when sleptons are produced in the decay 
of directly produced charginos-neutralinos. 

• Indirect slepton production from gluinos and squarks through charginos, neutrali- 

nos (domain 11), leads to a two {or more) leptons + E'™''* + jets signature. This 
production has the largest cross-section because of the strongly interacting g and g, 
and again sleptons can be singly produced. 

The largest contribution to the two leptons with no jets final state in charged slepton 
pair production comes from their direct decays to the LSP, II ^ I Xi- However, with 
decreasing mi/2 along an isomass curve the fraction of slepton direct decays to LSP 
decreases as shown in Fig. 2.9a. Figure 2.9b The a(/L/L) xB(2/ + invisible) for directly 
produced left sleptons and Fig. 2.9c shows the same for right sleptons: a{lRlji)x'B{2l + 
invisible). As right sleptons decay directly into LSP's the cross-section times branching 
ratio is constant along an isomass curve. One can see that the topological cross-section 
xB(2/ -|- invisible) is larger than (t(/l/l) xB(2/ -|- invisible) in almost the entire 
(mo,mi/2) plane. 

For sneutrino pair production, the main contribution to the 2 leptons and no jets 
topology comes from one sneutrino decaying directly to the LSP and the second to xf 
or X2 followed by their leptonic decays. The sneutrino decay branching ratios are shown 
in Figs. 2.10a-c. The o"(z/z/) xB(2/ -|- invisible), shown in Fig. 2.10d, is limited along the 
mi/2 axis by the xf , X2 leptonic decay branching ratios (see Figs. 2.1 and 2.2), and along 
the thq axis by the sneutrino production cross-section. 

Figure 2.11 shows the total cross-section times branching ratio for two lepton+ EJp'^^^ + 
no jets event topology arising from all possible combinations of directly produced slepton 
pairs, IlIl, IrIr, i^i^- It drops rapidly with increasing rriQ and mi/2, from a maximum 
value of the order of a few pb down to a few fb for mo ~ mi/2 ~ 500 GeV. In almost the 
entire (mo, mi/2) parameter plane slepton pair production is dominated by right sleptons. 

To illustrate the additional contribution from xf xf direct production and the poten- 
tial contribution from xf X2 ^'^ ^'^'^ same-flavor opposite-sign leptons+ i?™*** + no 
jets event topology, we show o-(xf xf ) xB(xf xf ^ + invisible), o-(xf X2) xB(X2 ^ 
/f + /^) and cr(xf X2) xB(x2 + ^^) in Figs. 2.12a, b and c, respectively. One can see 
that the contribution from electroweak indirectly produced sleptons is not negligible, and 
in some areas of the (mo, mi/2) parameters space is even higher than for direct production. 
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Indirect slepton production from strongly produced g, q is larger still, when allowed, but 
is (almost) always accompanied by jets in the final state, so will not contribute unless the 
jets somehow escape detection. 

2.3 Gluinos and squarks 



In mSUGRA the gluino mass is determined mainly by mi/2 [0: 

Mg ^ 2.5mi/2 (8) 

The five squark fiavors {d, u, s, c, b) with their left and right chiral states are assumed 
to be mass degenerate, giving altogether ten degenerate squark fiavors with: 



^ ^ml + 6m2/2 (9) 

The II and tn are however treated differently. Through their mixing they give rise to two 
mass eigenstates ti and ^2 with ti most often significantly lighter than all other q states, 
which requires a separate discussion in terms of both, production and decay. Mixing 
also introduces a small non-degeneracy among the hi and 62 states, but this is usually 
neglected compared to behavior of the stop system. 

At LHC energies, the total SUSY particle production cross-section is dominated by 
strongly interacting gluinos and squarks, which through their cascade decays, (Figs. 2.13 
and 2.14) can produce many jets and leptons, with missing energy due to at least 2 escap- 
ing LSPs and possibly neutrinos. Figure 2.15 shows the total squark/gluino production 
cross-section contours versus mo and mi/2- Figure 2.16 shows the (mo,mi/2) parameter 
space subdivided into several domains corresponding to the characteristic kinematic con- 
figurations of gluino and squark events and the predominant sources of lepton production. 
Isomass contours for squarks, gluinos, light and pseudoscalar Higgses are also shown in 
Fig. 2.16 for tan/3 = 2, = 0, /i < 0, with h masses calculated at the two-loop radiative 
correction level. 

The first domain is characterised by nig > rriq and m^o > mj. This means that gluinos 
decay to quark-squark pairs but the contribution from gluino production to the total 
squark/gluino cross-section is small compared to the squarks. Both the X2 xf this 
domain decay via sleptons (see Figs. 2.2 and 2.3) giving rise to a significant number of 
leptons in the final state. This region extends to higher values of mo at fixed mi/2 than 
just two-body decays of X2 xt to sleptons, since the three-body decays (e.g. X2 ~^ 
Xi l^l^) still remain even when two-body modes X2 ^ Xi^ and — > XiW^ are open. 

In the second domain of the (mo,mi/2) plane, the gluino mass is still somewhat larger 
than the squark masses, so that gluinos still decay into quark-squark pairs. The yield of 
stop-top final states is however higher than in the first domain since the contribution of the 
gg process in the total g/q production cross-section is larger and the g ^ it branching 
ratio is also higher. This results in an increase of the jet multiplicity and in a higher 
fraction of &-jets in the final state. The X2 ^^"^ Xi are now lighter than the sleptons, and 
they decay to LSP -|- h/W^, respectively, which means a reduced yield of isolated leptons 
when compared to the first domain. 

In the third characteristic region of the (mo,mi/2) plane, the gluino mass is smaller 
than the mass of squarks, except for ti and bi, and hence the stop-top final state dominates 
in overall production. The cross section for qq is significantly smaller than that for gg. 
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All this gives rise to a significant jet multiplicity in the event and some increase of the 
yield of isolated leptons, as a result of t ^ Wb Ivh decays, and an enrichment in the 
fraction of 6-jets. Compared with the previous domain, nothing changes from the point 
of view of the X2 ^-J^d Xi decays. 

The fourth domain is characterised by very massive squarks compared to gluinos, thus 
gg production dominates and the gluino decays predominantly via three-body final states 
(mainly LSP + tt). This increases further the jet (and 6-jet) multiplicity of the event. 
Again the decays of the X2 ^-^id Xi ^-^^ the same as in domains 2 and 3. 

There is a fifth domain which has approximately the shape of a band along the mo 
axis, under domains 2-4, for mi/2 < 180 - 190 GeV. In this region the masses of xf and 
X2 a-re smaller than Xi + W and x? + h, respectively, and smaller than the slepton and 

sneutrino masses, thus 3-body decays of X2 XiQ^/Xi^^l ) xf XiQl'/Xi^^^) 
take place. However, the yield of leptons does not differ significantly from that of domain 
1. 
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mSUGRA parameters: tan p= 2, Aq= 0, |j, < 
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Figure 2.1: Isomass contours for: a) light {xi, Xi, X2) and b) heavy {xt, Xs^ xl) 
charginos / neutrahnos. 
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mSUGRA parameters: tanp = 2, A^. = 0, ji< 




Figure 2.2: Chargino decay branching ratio versus (mo,mi/2): a) xt ~^Xi ^) b) xt 
^It^^ X? ^, c) xf -^Ir^^ X? ^ and d) xf ^ X? ^ xS ^■ 
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mSUGRA parameters: tanp = 2, A^. = 0, |i< 
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mSUGRA parameters: tanP = 2, Aj. = 0, |i< 




mSUGKA parameters: tanp = 2, Aj. = U, |i< U 





Figure 2.5: X2 inclusive production and decay: a) branching ratio B{g X2+ x); 
b) B{ul —>■ X2+ c) B(x2 — + invisible) and d) X2 inclusive production 
cross-section times branching ratio into 
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mSUGKA parameters: tanp = 2, Aj. = U, |i< U 




Figure 2.6: xf inclusive production and decay: a) branching ratio B(^ xf+ x), 
b) B{ul ^ xf+ X), c) B(x^ invisible) and d) xf inclusive production 

cross-section times branching ratio into Z^. 
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tanjS=10, /i<0, Ae=0 tan^=)0, iu.>0, Ao=0 




Figure 2.7: The sum of X2 ~^ Xil~^l a^id X2 ^ ^ Xi^^^ branching ratios for 
various values of tan/3 = 10 and 35 for both signs of /i. 
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Slepton isomass contours: 




Figure 2.8a: Isomass contours for a) left charged sleptons, b) right charged sleptons 
and c) sneutrinos in mSUGRA parameter space (mo,mi/2) for tan/? = 2, Aq = 0, /j, < 0. 
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Figure 2.8b: Slepton isomass contours for m^-^ = mj^ = m^o=150 GeV corresponding 
to lines 1-3 respectively and mj^ = mj^ = m^o=400 GeV, lines 4-6 respectively. Line 
7 separates domains I and II of {mo, 1711/2) parameter space; in domain II charginos, 
neutr alines can decay to left sleptons. 



17 





Figure 2.9: a) Branching ratio of /l — ^ Xi + b) cross-section times branching ratio 
for 2leptons + i^™*'' + no jets events from direct left slepton pair production; c) cross- 
section times branching ratio for 2leptons + E'j?*** -|- no jets events from direct right 
slepton pair production. 
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Figure 2.10: Decay branching ratios of sneutrinos: a) — >■ x5 + ^) i> — >■ + ^ 
c) — > + as a function of (mo, rtii/2)] d) cross-section times branching ration for 
2leptons + E^^^^ + no jets events topology from sneutrino pair production. 
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Figure 2.11: Total cross-section times branching ratio for 2 leptons + i?™"* events 
arising from all possible combinations of directly produced slepton pairs, namely Il^l, 
IrIr and i>i>. 
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Figure 2.12: Cross-section times branching ratio for a) xtxt direct production fol- 
lowed by leptonic decays of at least one of the charginos; b) single sleptons production 
from directly produced xtx^ with decays to left sleptons X2 ^ ^lI (b); c) single sleptons 
production from directly produced xfx2 with X2 decays to right sleptons X2 IrI (c). 
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Figure 2.13: Typical decay modes for the massive (1 TeV) gluino illustrating the 
variety of possible cascade decays. 
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Figure 2.14: Typical decay modes of the massive (0.95 TeV) squark. 
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Figure 2.15 Gluino/squark production cross-section versus (mo,mi/2)- 
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Figure 2.16: Domains of the (mo,mi/2) parameter space with charachteristic predom- 
inant decay modes. Isomass contours for squarks, gluinos, hght and pseudoscalar higgses 
are also shown. The shaded region near the mi/2 axis shows the theoretically forbidden 
region of parameter space, and a similar region along the mo axis corresponds to both, 
theoretically and experimentally excluded portions of parameter space. 
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Figure 2.17: X2 ~^ Xi^ branching ratio versus (mo, mi/2) for fixed Aq — and different 
values of tan/3 and /i. 
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3 Experimental signatures considered 



As discribed earlier, the highest cross-section for high mass i?-parity conserving SUSY 
at hadron coUiders is due to squarks and/or gluinos which decay through a number of 
steps to quarks, gluons, charginos, ncutrahnos, sleptons, W, Z, Higgses and ultimately 
a stable LSP (lightest neutralino in mSUGRA), which is weakly interacting and escapes 
detection. The final state has missing energy (2 LSP's + neutrinos), a number of jets, 
and a variable number of leptons, depending on the decay chain. Due to the escaping 
LSP's, which appear at the end of each sparticle decay chain, the masses of sparticles 
can not generally be reconstructed explicitely. However, the sparticle production and 
decay characteristics, discussed in the previous section, lead to a number of specific event 
topologies, which should allow the discovery of SUSY in general, and specifically the 
separation of certain SUSY sparticles and processes from SM and other SUSY processes. 
The LHC must therefore not only be able to discover SUSY, if it is reahzed at electroweak 
scale, and if it has not yet been discovered at the Tevatron (which covers a very promising 
energy/mass range), but the LHC experiments sholud also be able to disantagle the 
various SUSY production mechanisms and reconstruct most of the sparticle spectrum, 
determine or limit the models/scenarios, constrain model parameters, etc. Our present 
understanding shows that all this is possible at the LHC, but a lot of methodical studies 
are still required. Usually, one characterizes the general SUSY signal significance by an 
excess over the expected SM event rates, while for specific SUSY processes the background 
includes besides SM processes other SUSY reactions contributing to the same final state 
topology. 

The following final states have been investigated in substantial detail in CMS: 

1) leptons + jets + Eip'^^^ final states with a variable number of leptons and jets 
and possible requirements on 6-jets. These channels provide the maximum mass and 
parameter reach of SUSY through the production of strongly interacting particles, which 
in their cascade decays can give rise to many leptons and jets. 

2) inclusive 2 leptons + + {jets) and 3 leptons + (E"™*^) final states. In these 
channels SUSY reveals itself through X2 inclusive production with subsequent decay - 
directly or via a slepton - into two leptons and For some regions of parameter space, 
the two lepton invariant mass with its characteristic shape, a sharp edge at the end point 
of the spectrum, allows the determination of the sparticle masses, including the slepton 
mass, and model parameters. These channels could well be the first indication/signature 
for SUSY production at the LHC. 

3) exclusive 2 leptons + no jets + i?™®* final state, which is enriched in direct 
(Drell-Yan) slepton pair production as discussed in section 2.2. The main issue here is 
to understand and keep under control both the SM background (with many processes to 
be considered) and internal SUSY backgrounds. This channel could allow discovery and 
study of slepton production. 

4) exclusive 3 leptons + no jets + i?™** final state, which is the signature for direct 
xf X2 production in a Drell-Yan process. It is therefore theoretically most reliable (along 
with slepton production). For this and the slepton study, to suppress the SM and internal 
SUSY backgrounds a good response to jets down to E'^f* ~ 30 GeV, with a calorimetric 
coverage up to \ri\ ~ 4.5, is required. This exclusive channel could also play a central role 
in a precise determination of the x5 ^ the SUSY dark matter candidate. 

For the present study, production of SUSY processes has been simulated using the 
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IS A JET 7.14 Monte Carlo event generator fT^j and the SM backgrounds by PYTHIA 

Before presenting the physics 



5.7 



19], both with the CTEQ2L structure functions [2C]. 
expectations, we discuss the CMS detector optimisation studies done specifically with 
SUSY searches in mind. 



4 Detector issues 

4.1 CMS detector optimisation for SUSY studies 

The design goals of CMS are to measure muons, electrons and photons with a resolution 
of < 1% over a large momentum range (<100 GeV), to measure jets with a resolution 
of 10% at Et = 100 GeV and to be highly hermetic, with a missing Et performance as 
required for SUSY searches. The central element of CMS is a 13 m long, 6 m diameter 
solenoid generating a uniform magnetic field of 4 T pT| , p2| . Fig. 4.1. The magnetic flux 
is returned through a 1.8 m thick saturated iron yoke instrumented with muon chambers. 
Muons are precisely measured in the inner tracker and are identified and measured in 
four muon stations inserted in the return yoke. Precision tracking in the muon stations 
is carried out with drift tube planes in the barrel and with cathode strip chambers in 
the endcaps The goal is to achieve a spatial resolution of ~100 fim and an angular 



accuracy on a local muon track segment of ~1 mrad per station. Excellent time resolution 
is needed to identify the bunch crossing (with a periodicity of 25 nsec). Muon stations 
therefore also include resistive plate chamber triggering planes with time resolution ~2 
nsec 

The inner tracking system of CMS is designed to reconstruct high-p^ muons, isolated 
electrons and hadrons over \r]\ < 2.5 with a momentum resolution of Apx/pr — 0.15pT © 
0.5% {pt in TeV). Silicon and gas microstrip detectors are used to provide the required 
precision and granularity. In the present design there are about 9x10^ MSGC detector 
channels, about 4x10^ Si microstrip channels and about 8x10^ Si pixel channels p| . 



For high momentum muons the combination of tracker and muon chamber measurements 
greatly improves the resolution: Apx/pr — 0.06 for a p ~ 1 TeV muon in \ri\ < 1.6 [ pT| . 

The calorimeter system of CMS is made of a high resolution lead-tungstate (PbW04) 
crystal electromagnetic calorimeter and a copper-scintillator hadron calorimeter behind 
it. The primary function of the electromagnetic calorimeter is to precisely measure elec- 
trons and photons, and in conjunction with the hadron calorimeter to measure also jets. 
The arrangement of crystals is shown in Fig. 4.1. In the barrel the crystals are 25 ra- 
diation lengths [Xq) deep and the lateral granularity is ~ 2 cm x 2 cm corresponding 
to At] X A(j) = 0.014 X 0.014. They are read out with Si avalanche photodiodes p5| . 
The electromagnetic calorimetry extends over \ri\ < 3. The total number of crystals is 
~ 1.1 X 10^. With a prototype PbW04 calorimeter in a test beam, an energy resolution 



oi aE/E ~ 0.6% has been obtained for electrons of E = 120 GeV p6[. Hadron calorime- 
try with large geometrical coverage for measurement of multi-jet final states and missing 
transverse energy (E^^^^) is essential in all sparticle searches, as it is the EJp^'^'^ which 
provides evidence for the escaping LSP's (lightest neutralinos) . The hadron calorimeter 
of CMS is made of copper absorber plates interleaved with scintillator tiles read out with 



embedded wavelength shifting fibers The readout in the 4 Tesla field is done with 



hybrid photodetectors ES]. The tiles are organized in towers (Fig. 4.1) giving a lateral 
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segmentation of Ar/ x A</) ~ 0.09 x 0.09. The hadronic resolution obtained in a test beam 
is (Je/E ~ 100%/-\/E' © 5% for the combined PbW04 and hadronic calorimeter system 
p9| . This central hadron calorimetry extends up to |?7| = 3.0. It is complemented in 
the forward region 3.0 < \ri\ < 5.0 by quartz-fiber "very forward calorimeters" (Fig. 4.1) 



27| , pQ] . Their function is to ensure detector hermeticity for good missing transverse 



energy resolution, and to extend the forward jet detection and jet vetoing capability of 
CMS which is essential in slepton, chargino, neutralino searches as discussed in the fol- 
lowing. Detector hermeticity is particularly important for processes where the physical 
(real) missing Et is on the order of few tens of GeV as is the case in h, H, A ^ rr, W 
^ lu, t ^ lub, t H^b Tub, etc. and in particular in slepton, chargino and neutralino 
searches connecting the LEP2/Fermilab and the LHC search ranges. 



4.2 HCAL optimisation and tail catcher 

A number of steps have been undertaken in the design of CMS to optimise its E^^^^ 
response in view of SUSY searches. These include the inclusion of very forward calorime- 
try, the addition of a tail catcher behind the coil |]HT]] and optimisation of the crack for 
the passage of services between the barrel and endcap calorimeters. We now briefly dis- 
cuss these issues. The central {\ri\ <3.0) calorimetry of CMS is complemented by very 
forward quartz-fiber calorimeters [0 covering the rapidity range 3 < |?7| < 5, Fig. 4.1. 
Figure 4.2 illustrates the effect of the forward energy flow containment and its effect on 
missing Et measurements. It shows the expected instrumental (fake) missing Et distri- 
butions for QCD di-jet events as a function of calorimetric coverage at large rapidities. 
Calorimetry extending up to jr^l ~ 5 reduces the fake (instrumental) EJp'^'^^ by an or- 
der of magnitude in the 20 - 120 GeV range. Furthermore, the missing Et resolution, 
(T(i?™**)/S£'T where T,Et is the calorimetric transverse energy sum, is also much reduced 
in the presence of very forward calorimetry, see Fig. 4.3. There is a second requirement on 
forward calorimetry which turns out to be essential if we hope to eventually extract slep- 
ton and chargino/neutralino signals. In the search for direct DY slepton pair production 
// XiXi, or for the associated direct (DY) chargino/neutralino X^X2~^ ^^'^Xi^^^ Xi 

production, leading to final states with two or three isolated leptons, no jets, and missing 
Et, it is essential to have the capability to recognise and veto on forward jets. This is 
needed to suppress the large backgrounds due to tt, gg, gq, gg and the associated pro- 
duction modes qx-, gx which would otherwise overwhelm the signals. The extension into 
the forward direction of the central jet vetoing capability provides the needed additional 
rejection factors to keep these backgrounds under control. 

Figure 4.4a shows for the case of direct X1X2 production, (discussed in more detail in 
section 8), the expected rejection factors against the main backgrounds as a function of the 
rapidity range and detection threshold over which the jets can be recognized. Figure 4.4b 
shows the rejection factors against internal SUSY and the SM tt backgrounds. Here the jet 
veto is applied after lepton isolation and an Et^^^ cut, thus the value of the rejection factor 
is much reduced. In both cases, the jet coverage needed to obtain sufficient background 
rejection is |?7| up to ^ 4 and the loss of signal acceptance due to the jet veto is typically 
~10%. 

The locations of the tail-catcher scintillator layers, two in the barrel region behind the 
coil and one in the endcap, are shown on the longitudinal cut through the CMS detector 
in Fig. 4.1. Figure 4.5 shows the depth in interaction lengths A of calorimetry in CMS, 
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and the total sampled depth including the tail catcher layers, as a function of rapidity ||32 |. 
Whilst the total calorimetric absorber thickness within the coil is somewhat marginal at 
rapidity 77 ~ for full hadronic shower containment, inclusion of tail-catcher layers allows 
hadron energy measurements with at least 10.5 interaction lengths everywhere. 

Figure 4.6 shows the effects on hadron energy measurements of the inclusion of the 
tail-catcher layers, specifically on the reduction of the low-energy tail in the response to 
(200 GeV) hadrons What is shown is the ratio of the (GEANT) simulated response to 
the incident energy for single pions in the rapidity range 0.33 to 0.86 with no tail catcher, 
with one, and with two tail catcher layers included in the readout. Mismeasurements 
of hadronic energies, and in particular the presence of low-energy tails in hadronic (jet) 
energy measurements is one of the main sources of fake instrumental missing Et- When 
due to insufficient shower containment in depth, inclusion of the tail-catcher layers cures 
this problem to a large extent. 

There are two other sources of hadronic mismeasurements and low-energy tails: dead 
areas and volumes due to detector cracks and matching the hadronic response of ECAL 
to HCAL. 

Special attention has been paid in CMS to the effects on hadronic energy measurements 
- and thus on missing Ex - of the main crack for services between the barrel and endcap 
calorimetry at rapidity ~ 1.2 - 1.5 (Fig. 4.1). The various configurations which have been 
studied for this crack are shown in Fig. 4.7 |T^. The aim was to minimize the degradation 
of the hadronic energy resolution and the level of the low energy tail for hadrons and 
jets straddling the region of the gap. The optimal choice, from the point of view of 
hadron response, weight of the cantilevered endcap hadron calorimeter, manoeuvrability 
of endcaps and margin of freedom to close the detector, is a conical endcap shape with 
a non-pointing crack at 52 degrees from the beam line labeled TP7 in Fig. 4.7. Figure 
4.8 shows the reconstructed energy response (GEANT) to single pions of 100 GeV when 
incident on and off the crack region, for three crack widths of 11, 18 and 25cm. With 
the present estimate of the volume of services a gap of 12 cm is needed, plus 2 cm for 
clearance. 

An essential ingredient to obtain optimal EJp^'^^ response of a detector is the linearity 
of its response as a function of incident hadron (jet) energy. The non-compensating 
mixed calorimeter of CMS - with a PbW04 crystal ECAL compartment followed by a 
Cu/scintillator sampling HCAL - has significantly different electromagnetic to hadronic 
(e/h) responses in the two parts, e/h ~ 1.8 in the ECAL vs. 1.2 in the HCAL. Linearity is 
more difficult to achieve in such a mixed calorimeter than in a more homogeneous system. 
To restore linearity as much as possible, and to compensate for the effects of dead material 
in the space between the ECAL and HCAL, weighting techniques for responses must be 
used between ECAL and HCAL, and within HCAL compartments and tail catcher layers. 
An essential limiting factor is the number of readout channels we can afford for the 
longitudinal HCAL tower segmentation. The optimal set of weighting factors must be 
determined on basis of test beam data and detailed simulations of responses to hadrons 



and jets, taking into account financial limitations. In the present HCAL design [2/ 



a good compromise between expected performance and cost is found by having a very 
shallow (1 scintillator plane) front HCAL tower readout and second deep (~ 5 A) HCAL 
segment readout. The first shallow HCAL segment just behind the ECAL detects hadron 
interactions in the ECAL - where their energy deposition is underestimated (e/h ~ 2) - 
and compensates by overweighting this layer. Present investigations show that with this 
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technique the expected non-hnearity should not exceed 4% for pions between 20 and 300 
GeV|2l. 



4.3 The role of the tracker in SUSY searches 

The inner tracking system of CMS is primarily designed to reconstruct high-p^ muons, 
isolated electrons and hadrons over \t]\ < 2.5 with a momentum resolution of Apt/pt ~ 
0.15pT©0.5% {pt in TeV). Hadrons must be reconstructed down to ~ 1 - 2 GeV, as lepton 
and photon isolation is a very important selection criterion for a number of physics signals, 
and in particular in the SUSY searches for sleptons, charginos and neutralinos. Another 
important task of the tracker is to measure track impact parameters allowing the detection 
and measurement of long-lived particles and the tagging of 6-jets. The main problem in 
tracking is that of pattern recognition. At a luminosity of 10^^ cm~^s~^, interesting 
events will be superimposed on a background of about 500 soft charged tracks within the 
rapidity range considered, coming from ~ 15 minimum bias events which occur in the 
same bunch crossing. To solve the pattern recognition problem, detectors with small cell 
sizes are required. Silicon and gas microstrip detectors provide the required precision and 
granularity to maintain cell occupancies below 1%, but the number of detector channels 
is large (~ 1.5 x 10''). The design of the CMS tracker is shown in Fig. 4.9. 

A track in the barrel part of the tracker, in its present design (summer 1997), first en- 
counters two layers of pixel detectors, then four layers of microstrip Si detectors of 67 and 
100 yum pitch followed by seven layers of 200 fim pitch microstrip gas chambers (MSGC; 
~50 yum precision) Alternate MSGC and Si microstrip layers are double-sided to 

allow determination of the z coordinate of a track by a small-angle stereo measurement. 
The inner cylindrical volume (r < 50 cm) with Si detectors will be kept at ~ —5° C 
temperature, whilst the outer (MSGC) part at ~ 18° C. Figure 4.10 shows the expected 
track momentum resolutions in the CMS tracker alone. As already mentioned, for high 
momentum muons the combination of tracker and muon chamber measurements greatly 
improves the resolution: Ap^/pT — 0.06 for a p ~ 1 TeV muon in \t]\ < 1.6 This 



requires that the relative alignment of the inner and outer systems be known to within 



~100 /im p3[. The pixel layers, at radii of 7.7 and 11 cm in the high luminosity pixel de- 
tector option, ensure precise impact parameter measurements, with an asymptotic (high 
momentum) accuracy of ajp = 23 fim in the transverse plane (Fig. 4.10), and of ~ 90 
/im along the z axis [^. For the initial "low luminosity" (< 10^^ cm~^s~^) running, the 
impact parameter performance is improved by having the pixel layers at 4 and 7.7 cm 
radii. 



4.4 Lepton isolation 

In the search for direct DY slepton pair production // XiXi or for the associated 

direct (DY) chargino/neutralino xtxt^ li^Xi^~^^~Xi production leading to a final state 
with two or three isolated leptons, no jets and missing Et, there are two essential instru- 
mental requirements allowing separation of the signal. The first is lepton isolation, to 
suppress the copious backgrounds from processes such as tt, Zbb, bb. The rejection factors 
expected on basis of tracker isolation criteria as a function of the px cut on the accompa- 
nying tracks are shown in Fig. 4.11 for low and high luminosity running conditions. The 
second requirement, as discussed previously, is on the detector capability to veto on jets. 
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4.5 Tagging of 6-jets in CMS 

As discussed in Section 2.3, an important ingredient for SUSY physics studies is the 
capabihty of the detector to tag 6-jets. Figures 2.13, 2.14 and 2.15 illustrate the many 
ways 6-jets can serve as final state signatures of 6, t, t production in the g/q cascades. 
Particularly important is the possibility to detect h — >■ 66 in ^/g — >■ Xi^ decay channel 
with S/B r-u \ and potentially with much lower integrated luminosity than required to 
observe h 77. Techniques for 6-tagging have been studied in detail in CMS. The 
inner part of the CMS tracker, in particular the two pixel layers, has been optimised 
in detector positioning and required resolution with 6-tagging performance as the main 



criterion . The expected 6-tagging efficiency and sample purity has been studied with 
several tagging algorithms based on track impact parameter measurements. Figure 4.12 
shows, for example, the expected tagging efficiency as a function of jet Et averaged over 
the pixel detector acceptance. For 6-jets in ti events it is ~ 30%, whilst the expected fake 
rate probability is ~ 1%. The tagging algorithm used for Fig. 4.12 requires > 3 tracks of 
Pt > 2 GeV with > 2 a significance for impact parameter measurements in the transverse 
plane. 

How much the expected tagging efficiency is dependent on the tagging algorithm 
and what is the possible trade-off between efficiency and sample purity is illustrated by 
Fig. 4.13 where several tagging algorithms have been exercised. Figure 4.14 shows how 
the tagging efficiency and purity depends on the radius of the first barrel pixel layer, 
changing it from 4 cm as possible for the low luminosity running to 7.7 cm as required 
to sustain the radiation damage at 10^^ cm~^s~^. There would be an important gain in 
tagging efficiency if the first pixel layer could be always kept at 4 cm radius, at a cost of 
replacing it (almost) every year. This is particularly important for channels like h bb 
which go with the square of the tagging efficiency. Whilst the tagging efficiency is a 
simple function of the pixel point resolution and the radial position of the barrel pixel 
layers, the mistagging probability is dependent on the non-Gaussian part of the measured 
impact parameter distribution and is more difficult to control as it depends on the overall 
pattern recognition performance and multiple scattering of the tracker. All the 6-tagging 
efficiencies discussed here rely on impact parameter measurements; further improvements 
in 6-tagging efficiency are possible, including explicit reconstruction of secondary vertices, 
lepton tags, etc. 

4.6 CMSJET — an approximate description of detector response 

Understanding of the detector response with an adequate simulation is a very important 
part of correctly estimating the SUSY discovery potential of the CMS experiment. We 
have used the CMSJET program, which is a fast, non-GEANT simulation package for 



the CMS detector [^, CMSJET provides a satisfactory approximate description of 
detector components and of the response to hadrons, leptons and gammas. It is a good 
compromise between performance speed and precision. The following aspects of CMSJET 
are relevant to this study: 

• Charged particles are tracked in a 4 T magnetic field. A reconstruction efficiency of 
90% per charged track with > I GeV and \ri\ < 2.5 is assumed. 

• The geometrical acceptances for fi and e are \ri\ < 2.4 and 2.5, respectively. The 
lepton momentum is smeared according to parameterizations obtained from full GEANT 
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simulations. For a 10 GeV Icpton the momentum resolution IS^tIvt is better than 1% 
over the full 77 coverage. For a 100 GeV lepton the resolution becomes ~ 1^5% depending 
on 77. We have assumed a 90% triggering plus reconstruction efficiency per lepton within 
the geometrical acceptance of the CMS detector. This value is probably pessimistic for 
muon tracks, but more realistic for electrons. 

• The electromagnetic calorimeter of CMS extends up to \!]\ = 2.61. There is a pointing 
crack in the ECAL barrel/endcap transition region between \ri\ = 1.479— 1.566 (6 ECAL 
crystals) with significantly degraded resolution. The hadronic calorimeter covers |?7| < 3. 
The Very Forward calorimeter extends from \ri\ = 3 to \ri\ = 5. The full granularity of 
the calorimetric system is implemented in CMSJET. The granularity, assumed energy 
resolutions and electronic noise of calorimeters are listed in Table 4.1. Noise terms have 
been simulated with Gaussian distributions and zero suppression cuts have been applied. 



Table 4.1: CMS calorimeter description 





T] region 


Arj X A(j) 


AE/E 


Noise (MeV) 


ECAL 


|r/| < 1.479 
1.566<|77|<2.00 
2.00<|7?|<2.35 
2.35< 7? <2.61 


0.0145x0.0145 

0.0217x0.0218 
0.0292x0.0291 
0.0433x0.0436 


5%/^/E®0.5% 


25, 3(7 zero sup. 


Crack 


1.479<|r?|<1.566 


0.0870x6.283 


50%/VE®2% 




HCAL 


|r7| < 2.262 
2.262<|7?|<2.61 
2.61<|??|<3.0 


0.087x0.087 

0.174x0.175 
0.195x0.349 


77 dependent 

parameterization; 
82%//Ee6.5% at 7?=0 


250, 2(7 zero sup. 


VF 


3 < |?7| < 5 


0.17x0.17 ^ 
0.175x0.175 


172%/ Ve e 9% 
for hadrons 


250, 2a zero sup. 



• e/7 and hadron shower developments are taken into account by parameterization 
of the lateral and longitudinal profiles of the showers. The starting point of a shower is 
fluctuated according to an exponential low. 

• Jets are reconstructed using a cone algorithm, with a cone radius R = 0.4 0.9 and 
variable transverse energy threshold on jets depending on event kinematics. 

• For the high luminosity (L = 10^^ cm~^s~^) study, event pile-up is taken into account 
by superposition of "hard" pile-up events on top of signal and background. These are 
represented by PYTHIA QCD jet events with pt > 5 GeV. The number of superimposed 
pile-up events has been fluctuated with a Poisson distribution having a mean value of 15. 

The CMSJET model of the CMS detector is depicted in Fig. 4.15. 
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Figure 4.2: Missing Et distributions for QCD di-jet events as a function of calorimetric 
coverage at large rapidities. 
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mo = 1 00 GeV, = 1 00 GeV 
m(g) = 294 GeV, m(q) = 270 GeV 




Figure 4.4a: Rejection factors which can be achieved in chargino/neutrahno search 
using a jet veto, as a function of jet acceptance and jet detection threshold. 
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Figure 4.4b: Jet veto rejection factors against squark/gluino and tt backgrounds in 
slepton searches in two-lepton final states. 
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Figure 4.5: Depth of calorimetry and total sampled depth including the tail catcher 
layers as a function of rapidity in CMS. 
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Effect of tail catcher layers on reduction of tail 
in response to pions 
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Figure 4.6: Effects on hadron energy measurements of tail catcher layers in the central 
and next-to-central CMS wheels i.e. in two adjacent rapidity ranges. 
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HCAL Barrel / endcap designs 
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Effects of crack size on response 
to pions in "conical endcap" design 
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Figure 4.8: Effects on hadron energy measurements of tlie crack for services between 
tlie barrel and tlie endcap calorimetry for tliree crack sizes. 
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Figure 4.9: Layout of the CMS tracker. 



41 




Impact parameter resolution of CMS 
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Figure 4.10: Momentum and impact parameter resolutions as a function of pt at 
various rapidities. 
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Figure 4.11: The rejection factors against tt, expected on, the basis of lepton tracker 
isolation criteria as a function of the cut on the accompanying tracks in a cone of 
i? = 0.3. 
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b - tagging in CIVIS 
Low luminosity " pixel configuration 
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Figure 4.12: Tagging efficiency for b- and c-jets and mistagging probability for light 
quark and gluon jets. The tagging algorithm requires > 3 tracks of > 2 GeV with > 
2a significance for impact parameter measurements in the transverse plane. 
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Figure 4.13: Tagging efficiency for 6-jets and mistagging probability for light quark 
and gluon jets for various tagging algorithms. 
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CMSJET 3.5 calorimetry model 
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Figure 4.15: The CMS calorimetry model in CMSJET 3.5 
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5 Leptons + jets+ Ej^^^^ channel — 
search for gluinos and squarks 

To investigate the maximum reach in parameter space for SUSY via the production of 
strongly interacting sparticles we study multilepton+multijet + EJp^^'^ final states, result- 
ing from the cascade decays of these sparticles . As mentioned previously, a complete 
mass reconstruction of gluinos and squarks is impossible due to the presence of escaping 
particles, two LSP's plus possibly neutrinos. So SUSY signal observability is based on 
an excess of events of a given topology over known (expected) backgrounds. Among the 
leptonic final states studied, including always at least 2 jets + i?™'^'' in the final state, 
are one lepton (1/), two leptons of the same sign (2/ SS), two leptons of opposite sign (2/ 
OS), and 3,4,5 leptons regardless the sign (3/, 41, 51). As discussed in the following, these 
topologies have been investigated both with and without explicit requirements on muon 
isolation. The jets + EJp^^'^ topology, with no isolated leptons, is not discussed, since it is 
the most sensitive one to detailed detector simulation, and preliminary studies show that 
the mass reach is inferior to one of the channels discussed here 0. 

In order to establish the limits of the SUSY reach in the (mo,mi/2) parameter space 
for various integrated luminosities, the signal was generated at more than 100 points. At 
each point the signal is generated as a mixture of all possible combinations of the strongly 
interacting sparticle production processes (gg, gq, qq) with masses and decay branching 
ratios corresponding to this point. The set of mo,mi/2 points where the signal samples 
were generated for tan/? = 2, Aq = 0, /i < is shown in Fig. 5.1. Figures 5.2a-d give, at 
four representative points in the parameter space, the main sparticle masses and decay 
chains with corresponding branching ratios. 

The backgrounds considered in this study are: tt, W + jets, Z + jets, WW, ZZ, 
ZW, Zbb and QCD (2^2 processes including bb). There could be some double count- 
ing between Z + jets and 7M) processes, but in practice this means more conservative 
background estimates. 

No pseudorapidity cuts are applied during background generation to avoid distortion 
of the event kinematics. All processes were generated in the following pT intervals: 40 - 
100, 100 - 200, 200 - 400, 400 - 800, > 800 GeV. This allows the accumulation of significant 
statistics in the high-p^^ range in a more economical way than without subdivision of the 
generation pT range . 

The kinematics of signal events is usually harder than for the SM backgrounds, in 
particular for the interesting regions corresponding to the maximum reach in terms of 
squark/gluino masses. The cross-sections of the background processes can however be 
higher by orders of magnitude, and the high-pr tails of different backgrounds can have a 
kinematics similar to that of the signal. Figures 5.3-5.5 illustrate kinematical distributions 
for the signal at three different points in parameter space and compare it to background. 
More specifically in Fig. 5.3 one can see the comparison between kinematical distributions 
for the signal point (mo = 200 GeV, mi/2 = 760 GeV) and the distributions for the 
sum of backgrounds in 21 SS, 21 OS and 3/ final states. Both signal and background 
histograms contain only events satisfying first level selection criteria described below. 
Only the hardest jet and lepton in the event are shown in distributions. The last histogram 
in each Figure shows the distribution of E^"^, the scalar sum of all relevant transverse 
energies for leptons, jets and E^^''^. It reflects the overall transverse energy flow in the 
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event, the event "hardness". This variable is also rather insensitive to the transverse 
energy flow induced by event pile-up since it accounts only for those calorimeter cells 
which are attached to a jet, and is thus a useful quantity to cut on in the high luminosity 
regime. 

The topology of signal and background events is already rather similar after the first 
level selection cuts, in terms of angular distributions and circularity, as can be seen in 
Fig. 5.3a, thus it is not very useful to apply explicit cuts on these variables. The difference 
between signal and background in the lepton px distributions is also not as pronounced as 
might have been expected, as signal leptons are produced in cascade decays, thus losing 
"memory" of the hardness of the original process. For extremely high masses of squarks 
or gluinos in the ~ 2 TeV mass range, the difference in the angular and pip distributions 
between signal and the total SM background is still significant and helps in obtaining an 
effective data selection. Ultimately, all the cuts have to be justified from the point of view 
of the best observabihty of the signal over expected background. 

Since the kinematics of the signal vary significantly with the variation of mo and mi/2 
parameters. Figs. 5.3-5.5, the set of different cuts is applied and optimised independently 
in each point investigated. These cuts arc listed in Table 5.1. The optimisation procedure 
takes into account the number of jets per event, and parameters of lepton isolation, as 
discussed in details in the next subsection. The Eip'^^^ cut can be incremented beyond 
the values hsted in Table 5.1, keeping all other cuts untouched, to take advantage of the 
large values this variable can take when the signal is close to the boundaries of explorable 
parameter space where masses of squarks and gluino are in the 1-2 TeV range. 

We find that requiring more than 2 jets (up to > 4-5 for 1, 10 fb~^, up to > 6 for 
100 fb~^) in the event is useful in many cases, especially in domains 3 and 4 as defined 
in Fig. 2.16. The exception is domain 1, where requiring > 2 jets (as well as requiring 
muon isolation for all cases except for 4 and 5 lepton final states) is sufficient to get the 
best signal observability. 



Table 5.1: The various cuts apphed to optimise the reach. 



CUT 


p^ (GeV) 


p^ (GeV) 


p^/"'''' (GeV) 


£;^^^*(GeV) 


E'^'^ (GeV) 


CI 


20 


10 


40 


100 




C2 


20 


10 


150,80,40 


150 


500 


C3 


40,20,20 


150,80,40 


150 


500 


C4 


40,20,20 


200,120,40 


200 


700 


C5 


50,30,20 


250,150,40 


200 


1000 


C6 


80,50,20 


300,200,40 


250 


1200 


C7 


100,70,20 


300,200,40 


250 


1200 



It is also found that some moderate cuts: 



• Circularity > 0.1 

• 5ip (hardest lepton, Eip^'') > 10-15 deg 

are useful to improve signal observability everywhere in 11 final states, and in 21 OS final 
states, primarily in the region mi/2 > 600 GeV for an integrated luminosity > 10 fb""^. 
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From the chosen multiparametric space of cuts and constraints, we look for those 
which provide the highest value of mi/2 with a 5a significance excess for signal at each 
probed ttiq point and for each lepton multiplicity final state investigated 0. 

An important issue in this study is the lepton isolation criterion. The term "isolated 
lepton" means satisfying simultaneously the following two requirements: 

• no charged particle with px > 2 GeV in a cone with defined radius in t] - (p space 
around the lepton, and 

• T^E^^" in a cone ring 0.1 < R < 0.3 around the lepton impact point has to be less 
than a specified percentage of the lepton transverse energy. 

Figure 5.6 shows the acceptance for the signal in two points of parameter space, 
"SIGNAL 1" at mo = 200 GeV, mi/2 = 700 GeV chosen in the first characteristic 
domain in Fig. 2.16; "SIGNAL 2" at mo = 1000 GeV, mi/2 = 600 GeV in the third 
domain, and the acceptance for background. The comparison is done for the three final 
state topologies, namely, two lepton of same sign (2/ SS), two leptons of opposite sign (2/ 
OS) and three leptons (3/). This is done as a function of the tracker isolation cone size 
R = \/6ri'^ + 6ip'^ and for the set of cuts denoted as C2 in Table 5.1. Both electrons and 
muons are required to be isolated. Each curve corresponds to a specified calorimetric cut: 
no cut (circles), E!j?'-° in the cone ring 0.1 < i? < 0.3 not exceeding 10% (squares) and 5% 
(triangles) of the lepton Ex- 
it is rather difficult to deduce from these curves what is in general the most advanta- 
geous isolation criterion. Tight isolation is not very useful in case of "SIGNAL 2" in 21 
OS and 3Z final states, where it leads to a significant loss of signal. At the same time, the 
background is not suppressed sufficiently, especially in case of 21 OS, where it is mainly 
due to ti production giving well isolated leptons through the decay chain t Wb — > lub. 
So, the isolation criterion could also in principle be optimised at every mo, mi/2 point and 
for each leptonic topology so as to obtain the best signal observability. In practice, since 
the identification of non-isolated electrons is experimentally difficult, a minimal isolation 
is always required, e.g. charged particle isolation in a cone of R = 0.1 and a requirement 
of E!j?'-° in R = 0.3 not exceeding 10% of the electron Et- We find that a still tighter 
isolation does not further improve the signal observability at all the mo, mi/2 points in- 
vestigated. As for the muons, we find that a gain of ~ 50-150 GeV in reach in mi/2 is 
possible if no isolation is applied at all on them. This is valid for all lepton topologies and 
almost everywhere in the mo,mi/2 plane, except in the first domain in Fig. 2.16. This is 
due to the copious production of 6-quarks in the signal, especially in domains 2-4, giving 
rise to a significant yield of non-isolated leptons from 5-decays. So except when explicitly 
stated, we do not apply in the following any isolation requirement on muons. 

The case of high instantaneous luminosity running requires a separate discussion of 
lepton isolation due to event pile-up deteriorating lepton (electron) isolation. However, as 
we use only a loose isolation criterion to begin with, the criterion can easily be changed 
with only minor loss of efficiency. Furthermore, with a good tracker it should be possible 
to separate most of the tracks originating from vertices other than the leptonic one. Even 
in case of confusion with tracks from event pile-up vertices, the probability of the signal 
lepton to be "moderately" isolated (as we have defined above) in spite of pile-up is still 
97 - 98% when averaged over the entire ECAL, as we calculated using PYTHIA pile-up 
(MSEL = 2) with an average number of pile- up events per bunch-crossing of n = 25. So 
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event pile-up at high instantaneous luminosity does not affect significantly the results of 
our study. 

Figures 5.7-5.9 show the final results of our study of squark and gluino observability 
and mass search in mo, mi/2 parameter space in the various multilepton + multijet + 
j^rniss gj-^g^i states for our standard set of parameters tan/3 = 2, Aq = 0, /i < 0. The results 
are shown for increasing values of integrated luminosity from 1 to 100 fb~^. All the 5(j 
boundaries shown are obtained relaxing muon isolation whenever possible, but keeping 
always the minimal electron isolation criteria needed for unambiguous instrumental elec- 
tron identification. The maximal squark and gluino mass reach is obtained in 1 lepton 
final states. With 10^ pb~^ the maximal gluino mass reach varies from ^ 2.5 to 2.0 TeV, 
depending on mo, whilst the squark mass reach varies between 2.1 - 2.5 TeV. What is 
important is that already with 10^ pb~^ the < 1 TeV domain can be thoroughly explored, 
and according to the "no-fine tuning" arguments this mass domain is the most plausible 
one. 

In Fig. 5.8 one can see that the cosmologically interesting region |]3^ is entirely con- 
tained within the domain explorable already with 10 fb~^ in all topologies considered. 
The Figure shows the relic neutralino dark matter density contours corresponding to VlK^ 
values of 0.25, 0.5 and 1, where VL is the neutralino density relative to the critical density 
of the Universe and h is the Hubble constant scaling factor {h ~ 0.7). The prefered mixed 
dark matter scenario corresponds to values of between 0.15 and 0.4. 

Figure 5.10 illustrates the gain in mass reach which can be obtained when non-isolated 
muons are taken into account in the analysis. Figure 5.11 shows the expected signal-to- 
background ratio at the ba boundaries in each topology. It is always bigger than 1, 
thus we may expect that the maximum reach will not be very sensitive to a possible un- 
derestimation of the background cross-sections. To estimate the influence of systematic 
uncertainties on the background cross-sections on the ^, q observability contours, we in- 
creased all the background cross-sections by a factor of 2, and the W/Z + jets background 
by a factor of 5. Figure 5.12 shows how the 5cr significance contours are modified under 
these conditions for an integrated luminosity of 10 fb~^. The loss of reach is small in all 
cases. 

We also checked the sensitivity of our results to changes in some of our basic param- 
eters. The change of the sign of /x, keeping tan(3 = 2,Aq = 0, affects primarily charginos 
and neutralinos. It changes their masses, decay modes and branching ratios, and these 
alter the final state kinematic and topology. In the squark sector, a positive sign of 
decreases the ti mass by ^ 30 - 50 GeV compared to the /i < case. Consequently, final 
state leptons coming from stops and top quarks will be slightly "softer", thus requiring 
lower threshold cuts and causing modified background levels compared to the /i < case. 
Figure 5.13 gives the maximum expected reach of the CMS detector for squark and gluino 
searches in different multilepton final states for 10^ pb~^ with /i > j^. The reach is 
again greatest in the single lepton channel, and extends up to gluino and squark masses of 
~ 2 TeV. Relaxation of the muon isolation requirements increases the reach by 50 GeV 
in the single lepton case, and by ~ 200 GeV in the trilepton case. The 4Z channel with 



10^ pb~^ has an interesting explorable domain only if muons are not isolated ||38|. Figure 
5.14 gives the same results for Li^t = 10^ pb~^ [^]. With this luminosity and with full 
muon isolation (not shown), the 4/ maximal reach region in fact becomes "a band of 
observability" approximately between mi/2 ~ 200 GeV and mi/2 ~ 350 GeV pS. The 



lower boundary is mainly caused by a decrease of the branching ratio for g — > Xi't't- For 
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mo < 200 GeV the region of observability extends up to mi/2 ~ 650 GeV, as here there 
is a high yield of leptons from chargino and neutralino decays. 

The main conclusions of this particular study are the following: within the mSUGRA 
model, the SUSY signal due to production of q and g would be detectable as an excess 
of events over SM expectations up to masses rriq ~ rrig ~ 1 TeV with only 10'^ pb^^. The 
ultimate reach, for 10^ pb^^, would extend up to ~ 2500 GeV for small mo 400 GeV) 
and up to ~ 2000 GeV for any reasonable value of mg 2000 GeV); squark masses can 
be probed for values in excess of 2000 GeV. This means that the entire plausible domain 
of EW-SUSY parameter space (for tan/? = 2) can be probed with these lepton{s) + jets 
_l_ j^miss gj^g^i states. Furthermore, the S/ B ratios are > 1 everywhere in parameter space 
(with the appropriate cuts) thus allowing a study of the kinematics of g, g production and 



obtaining information on their masses The cosmologically interesting region Vth < 
1 can be entirely probed with an integrated luminosity not in excess of ~ 10'^ pb~^, thus 
with a very large safety margin. 
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Figure 5.1: The set of generated data points in the (mo,mi/2) mSUGRA plane for 
tan/3 = 2, Ao = 0, /i < 0. 
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Figure 5.2: Some typical points in parameter space. 
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Figure 5.3a: Comparison of the signal and background kinematics at one point in 
mSUGRA parameter space, after requiring two leptons of opposite sign (plus jet and 
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2 1 SS case, both leptons are isolated 
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Figure 5.3b: Comparison of the signal and background kinematics in the same point 
in in Fig. 5.3a, but for 21 SS final state topology. 
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Figure 5.3c: Comparison of the signal and background kinematics in the same point 
in in Fig. 5.3a, but for 3/ final state topology. 
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= 500 GeV , niy^ = 500 GeV 

2 1 SS case, both leptons are isolated 
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Figure 5.4: Comparison of the signal and background kinematics for a second 
point in parameter space. 
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Figure 5.5: Comparison of the signal and background kinematics for a third point 
in parameter space. 
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Figure 5.6: Effect of the lepton isolation criterion on signal and background, at two 
different points in parameter space: "SIGNAL 1" is ttiq = 200 GeV, mi/2 = 700 GcV, 
"SIGNAL 2" is mo = 1000 GeV, mi/2 = 600 GeV, as a function of Risoi, the tracker 
isolation cone size Risoi = \/Srf~\~S(p^ . The results are shown for three choices of calori- 
metric isolation parameter (sum of calorimeter energy deposit in a cone ring 0.1 < < 0.3 
around the lepton impact point) and for three event topologies. 
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Figure 5.7: 5cr gluino/squark contours in various lepton multiplicity final states for 1 
fb~^ and our standard set of model parameter values tan/5 = 2, Aq = 0, < 0. Squark 
and gluino isomass curves shown by dash-dotted lines. 
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Figure 5.8: 5cr contours as in Fig. 5.7, but for 10 {h~^. The relic neutralino dark 
matter density contours of Qh"^ = 0.25, 0.5 and 1.0 are also shown. The region of Qh"^ 
between 0.15 and 0.4 is prefered by mixed dark matter models |3^; the flh'^ = 1 contour 
is the upper limit. 
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Figure 5.9: Ultimate squark and gluino mass reach, in terms of a > 5o" observation 
contour for 100 fb~^, in various lepton multiplicity final states . 
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Figure 5.10: Influence of muon isolation on the 5a contours in various topologies for 
10 fb-i. 
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Figure 5.11: Signal to SM background ratio on the 5a reach boundaries in in various 
lepton multiphcity final state topologies, for 10 fb~^. 
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Figure 5.12: Sensitivity in the systematic uncertainties in background cross section on 
the 5cr contours in various lepton multiphcity final state topologies, for 10 fb~^. 
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Figure 5.13: Maximal squark/gluino reach for fi > for 10 fb~^ in different final state 
topologies. The relic neutralino dark matter density contours of Qh"^ = 0.25, 0.5 and 1.0 
are also shown for this set of model parameters. 
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6 Search for next-to-lightest neutrahno 



Leptonic decays of the second-lightest neutrahno, X2 have a useful kinematical feature: 
the dilepton invariant mass spectrum has a sharp edge near the kinematical upper limit 
with MJ+p = M^o — M^o in the case of direct three-body decays ^^^~Xi 
MJ^P = J{M% - M?){M? - M?o)/Mr in the case of two-body cascade decays x^-^ 

^^^L,R ~^ f'^^ Xi- This feature was first discussed in Ref. for xfx2 direct (EW) 
production. On the other hand, at LHC, the X2 is abundantly produced from gluinos and 
squarks over whole (mo,mi/2) plane, see e.g. Figs. 2.5a-b, and this source of X2 is much 
more prolific than direct EW production. Therefore, we propose to perform an inclusive 
search for X2 ^^'^ ^se the specific shape of the dilepton mass spectrum as evidence for 
SUSY HTI, PI . For the SM background suppression, besides two same-flavor opposite-sign 



leptons, one can ask for an additional signature characterizing SUSY. This can be: missing 
energy taken away by escaping LSPs, or additional jets coming from gluinos and squarks 
cascade decays, or an extra high-p^ lepton from charginos, neutralinos, sleptons, vector 
bosons, or 6-quarks copiously produced in SUSY. Depending on the sparticle masses and 
predominant decay modes one of these extra requirements can turn out to be more advan- 
tageous than the others, or to be complementary. In the following we discuss the inclusive 
2 leptons + E^^^^ -\- (jets) and 3 leptons -\- (i?™*^) channels and determine the region 
of {mQ,mi/2) parameter space where the dilepton invariant mass edge is visible in these 
final states. 



6.1 Inclusive 2 leptons + E'^*'^^ + (jets) channel 

The most straightforward signature for selecting the X2 decays is provided by the topology 
with two same-flavor opposite-sign leptons accompanied by large missing transverse energy 
and usually accompanied by a high multiplicity of jets. Here we concentrate on the 
2 leptons -\- i?™'^* -|- (jets) channel, where the flnal state leptons are electrons and muons 

EM- 

The largest SM background is due to tt production, with both W^s decaying into 
leptons, or one of the leptons coming from a W decay and the other from the 6-decay 
of the same t-quark. We also considered other SM backgrounds: W -\- jets, WW, WZ, 
Wtb, bb and rr-pair production, with decays into electrons and muons. Chargino pair 
production xtxi is the largest SUSY background but gives a small contribution compared 
to the signal. 

To observe an edge in the M;+;- distributions with the statistics provided by an inte- 
grated luminosity Lint = 10^ pb~^ in a signiflcant part of the (?77.o,wii/2) parameter plane, 
it is enough to require two hard isolated leptons (p^'^ > 15 GeV) accompanied by large 
missing energy, i?™*'' > 100 GeV. Our criterion for observing an edge in the Mi+i- distri- 
bution contains two requirements: {Nev — Nb)I\/Nev > 5 and {Nev — Nb)/Nb > 1.3, 
where Nev is the number of events with Mi+i- < MJ+p, and Nb is number of the ex- 
pected background events. Fig. 6.1a shows the invariant mass spectra of the two leptons 
at various (mo,mi/2) points for tan/5 = 2, fi < 0. The observability of the "edge" varies 
from 77a and signal to background ratio 31 at point (200,160) to 27a and a signal to 
background ratio 2.3 at point (60,230). The appearance of the edges in the distributions 
is sufficiently pronounced already with Lint = 10'^ pb"^ in a significant part of (mo,mi/2) 
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parameter plane. The edge position can be measured with a precision of ~ 0.5 GeV. 

With increasing mo and mi/2 the cross-sections are decreasing and therefore higher 
luminosity and harder cuts are needed. To achieve maximum reach in mi/2 with L^t = 
10"^ pb~^ for points where X2 has two-body decays via II {mQ < 0.45 ■ mi/2), a cut up to 
j^miss ^ ggg QgY jg ncccssary to sufficiently suppress the background. For points with 
large mo, in the region where X2 has three-body decay (mo ^ 0.5 ■ mi/2, mi/2 ^ 200 GeV) 
the transverse momentum pt of the leptons and E^^'^^ are not very large, but there are 
more hard jets due to gluino and squark decays. Thus for these points we keep the 
same cuts for leptons and missing energy as before [prj^'^ > 15 GeV, E^'^'^'^ > 100 GeV) 
and require in addition a jet multiplicity Njet > 3, with energy El'' > 100 GeV, in the 
rapidity range | rjjet \< 3.5. To optimise the edge visibility we also apply an azimuthal 
angle cut, A0(/"'"/~) < 120°. This jet multiplicity requirement is also helpful for the first 
domain, where X2 decays to leptons through (mo ^ 0.45 ■ mi/2). Right sleptons are too 
light to provide large lepton pt and i?™*^, and to use cuts on pip and E^^'''^ alone is not 
very advantageous. With Lint = 10^ pb~^, to suppress the background at larger accessible 
^0, ^1/2 values, we have to require at least 2 or 3 jets, depending on the mo, mi/2 region 
to be explored. Fig. 6.1b shows invariant dilepton mass distributions at some (mo,mi/2) 
points close to maximum reach with Lint = pb~^ and Lint = 10^ pb~^ respectively. 

Figure 6.2 shows the inclusive dilepton spectrum for larger values of tan/5 = 10 and 
35. Despite decreasing branching ratios as visible in Fig. 2.7, the dilepton (e"'"e~, n'^fJ^') 
edge is visible for low mi/2 values already with 10^ pb^^ |4^ . 

The regions of the (mo, rni/2) parameter plane where an edge in the Mi+i- spectra 
for tan/5 = 2, /i < can be observed at different luminosities are shown in Fig. 6.3. In 
Fig. 6.4 we show separately the three domains where an edge due to X2 ~^ ^^Xi? ^R^ and 
IlI decays can be observed at Lint = 10^ pb~^. One can notice a small overlapping region, 
where we expect to observe two edges, due to X2 ~^ Xi and to X2 ~^ ^r^^ ~^ Xi 
decays {case 1). With increasing luminosity and correspondingly higher statistics, this 
overlapping region increases, see Figs. 6.5 and 6.6. These plots show the same as Fig. 6.4, 
but for Lint = I0'^pb~^ and L^nt = 10^ pb~^, respectively. 



6.2 Inclusive 3 leptons channel 



In this channel an extra high-p^- lepton is required to suppress the SM backgrounds ^ . 
This can be either an isolated lepton, e.g. from the xf decay which has a similar behavior 
to X2 (Fig- 2.6), or a non-isolated lepton from the abundant production of high transverse 
momentum 6-jets in SUSY events. Thus in the inclusive 3 lepton channel we require: 

• two opposite-sign, same-fiavor leptons (/x or e) with pt > 10 GeV; 

• a third "tagging" lepton with pt > 15 GeV; 

• lepton isolation: if there is no track with Pt > 2 (1.5) GeV within AR = y/Ar]'^ + A(f)'^ < 
0.3 about the lepton direction, it is considered as isolated. 

The lepton px thresholds can be varied depending on the region of (mo, mi/2) studied 
and in some cases leptons are not required to be isolated. In the selected 3 lepton events 
we reconstruct the invariant mass Mi+;_ of a lepton pair with the same fiavor and opposite 
sign. When several /"'"/" combinations per event are possible, the one with minimal sepa- 
ration ARi+i- is chosen. Any significant deviation from the expected SM l^l~ invariant 
mass spectrum then provides evidence for SUSY. 

We stress that no other requirements, such as missing transverse energy, jet activity. 
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etc., are imposed in the inclusive 3 leptons channel. The price for this simplicity is a 
decrease of signal acceptance due to the requirement on the third lepton, and somewhat 
higher px of leptons needed to sufficiently suppress the SM backgrounds. Thus there 
is complementarity between this search which relies on lepton (e, /x) measurements and 
the tracker for lepton isolation and momentum measurement, and the search, in inclusive 
2 leptons + EJp^^^ channel which relies on lepton detection and overall calorimetry, HCAL 
in particular, for i?™'^^ measurement. 

The main SM sources of three leptons considered are production of WZ, ti, TIL and 
Z66. The expected SM background to 3 leptons final states at an integrated luminosity 
of Lint = 10^ pb~^ is shown in Fig. 6.7. Near the Z mass the main contribution comes 
from the WZ production, while tt production dominates at lower invariant masses. 

Dilepton spectra for the mSUGRA "Point 4" (see Table 1.1), superimposed on the SM 
background, are shown in Fig. 6.8. The number of events corresponds to an integrated 
luminosity of Lint = 20 pb~^, i.e. just a few hours of initial LHC running. The pt 
threshold on all three leptons is 15 GeV. In Fig. 6.8a no isolation requirement is applied, 
whereas in Fig. 6.8b the two leptons entering the dilepton mass distribution are isolated. 
The contribution from the SM background is negligible in the latter case. In both cases 
the specific shape of the /"*"/" distribution, with its sharp edge, reveals X2 production. 
At this mSUGRA point the X2 ^ three-body decay mode with a branching ratio 
B(X2 ^ + Xi) = 0.32 (/ = e, /x) and the masses are M^o = 97 GeV, M^o = 45 GeV. 
The sharp edge is situated at the expected value of = 52 GeV. 

Figure 6.9 shows the Mi^i^ distribution for mSUGRA "Point 1" (see Table 1.1) with 
an integrated luminosity of Li^t = 10^ pb~^. The pr thresholds for leptons are 15, 15, 30 
GeV and all three leptons are isolated. At this mSUGRA point the X2 possesses a two- 
body decay mode with the branching ratio B(x2 ^_r+0 — 0.24. The sparticle masses are 
M^o = 231 GeV, M^o = 122 GeV, Mr = 157 GeV. The "edge" of the distribution is sit- 
uated near 108 GeV, as expected. We define a signal significance as S" = Ns/y/Ns + Nb, 
where Ns and Nb are the numbers of signal and background events, respectively, in a 
mass window below the edge. In this particular case signal and background events are 
calculated in a mass window of 100 GeV< Mi+i~ < 110 GeV, Ns = 140 and A''^ = 39, 
resulting in signal significance of S = 10.5. Nb consists of 18, 8, 11 and 2 events coming 
from WZ, ti, TIL and Z66 production, respectively. 

An example for the mSUGRA "Point 5" is shown in Fig. 6.10. The cuts are the same 
as in the previous case. At this point X2 has three-body decay modes with B(x2 + 
Xi) = 0.06 and the masses are M^o = 124 GeV, M^o = 73 GeV. The observed value of 
M™"5 is close to the expected value of 51 GeV. In a mass window of 20 GeV< Mi+i- < 52 
GeV the numbers of signal and background events are Ns = 588 and Nb = 402 giving a 
signal significance of S" = 19. Nb consists of 113, 139, 65 and 85 events coming from WZ, 
ti, TIL and Z66 production, respectively. 

By the same procedure the mSUGRA parameter space (mo,mi/2) was scanned for 
fixed tan/3 = 2, = and /x < p. Figures 6.11a-c show the M;+;_ distributions for 
three different mSUGRA points with the optimal thresholds on lepton p-p and isolation 
requirement in inclusive 3 lepton events. Generally, with increasing mi/2 the observation 
of the edge becomes more difficult due to the rapidly decreasing gluino and squark produc- 
tion cross-sections. Nonetheless, this spectacular structure in the dilepton mass spectrum 
refiecting decays, may be among the first ways through which presence of SUSY may 
reveal itself. It may appear with very modest statistics, in a non-sophisticated analysis 
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and with only modest demands on detector performance. The precision for "edge" mea- 
surement in most of the region, where it is detectable is expected to be < 0.5 GeV. Note, 
that the Z peak seen in Figs. 6.8-6.11 serves as an overall calibration signal; it allows to 
control the mass scale as well as the production cross-section. 

Asking for a statistical significance S > 7 with at least 40 signal events calculated 
in a mass interval > 10 GeV below the edge the domain explorable with an integrated 
luminosity of Lint = 10^ pb~^ is shown in Fig. 6.12. It entirely covers the cosmologically 
prefered region 0.15 < Qh"^ < 0.4 P7|; here Q is the ratio of the relic particle density to 
critical density, h is the Hubble constant scaling factor and bounds are obtained assuming 
Xi is the cold dark matter particle in the mixed dark matter scenario for the Universe, 
varying h over the allowed range 0.5 < h < 0.8. 

For some regions of mSUGRA parameter space MJ+p is close to, or even hidden by 
the Z signal. At these points the accuracy of an edge measurement is ~10 GeV. Applying 
a cut on missing transverse energy (and/or jets) suppresses the contribution from the SM 
Z production, improving signal visibility. 



6.3 Inclusive 3 leptons + E"™*^ channel 

With increasing mi /2 , gluino and squark masses increase and the missing transverse energy 
becomes larger. Requiring £^™'^*> 200 (300) GeV rejects most of the SM background 
leaving a big fraction of signal events. Figures 6.11d-e show the dilepton spectra for 
various mSUGRA points in the 3 lepton + E"™**^ final states. In the region of high 
mi/2, dileptons are mainly produced in the cascade decays of X2- The third, "tagging" 
lepton predominantly comes from xf which also has cascade leptonic decays. Thus in this 
region of parameter space a large fraction of 3 lepton + E^^''^ events are in fact inclusively 
produced xfx2 pairs. Therefore all three leptons are required to be isolated. The lepton 
Pt thresholds are chosen to be quite asymmetric: 50, 25, 10 GeV, to account for the 
cascade nature of the X2 ^^cl xf decays. The region of mSUGRA parameter space, where 
the X2 is observable from the dilepton mass distribution shape in the 3 lepton + E'™'^* 
events at Lint = 10^ pb~^ luminosity is shown in Fig. 6.12 [^]. The reach extends up to 
mi/2 ~ 900 GeV and a detectable edge is seen as long as a-B > 10^^ pb (see Fig. 2.5). 

We conclude that 

• Observation of an edge in the dilepton invariant mass spectrum refiects the produc- 
tion of X2 hence establishes the existence of SUSY; this observation in some cases is 
possible with very small statistics and could be the first evidence for SUSY at LHC. 

• With no great demand on detector performance a large portion of mSUGRA parame- 
ter space including the cosmologically prefered domain can be explored with an integrated 
luminosity as low as Lint = 10'^ pb~^. 
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Figure 6.1a: Invariant mass distribution of two same-flavor, opposite-sign leptons at var- 
ious (mo,mi/2) points for tan/3 = 2, < with Lint = 10'^ pb^^ SM background is also 
shown (dashed hne). 
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Figure 6.1b: Invariant mass distribution of two same-flavor, opposite-sign leptons at 

(mo, ^1/2) points close to the experimental reach at corresponding luminosities Lint = 10^ 
pb~^ and 10^ pb~^ for tan/? — 2, /j, < 0. SM background is also shown (dashed line). 
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Figure 6.2: Invariant mass distribution of two same-flavor, opposite-sign leptons at various 
"^0^^1/2, tan/5 and /i points for Lint = 10'^ pb"^. SM background is also shown. 
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Figure 6.3: Observability of edges in invariant dilepton mass distribution with luminosities 
10^ pb~^ (dashed line), 10*^ pb"^ (solid line) and 10^ pb^^ (dashed-dotted line). Also shown 
are the explorable domain in sparticle searches at LEP2 (300 pb~^) and the Tevatron (1 
fb~^), theoretically and experimentally excluded regions. 
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Figure 6.4: Domains where the observed edge in the Mi+i- distribution is due to the 
decays X2 ~^ ^t^"^ ~^ Xi^'*'^" (sohd hne), X2 ~^ ^r^^ ~^ Xi^"*"^" (dashed-dotted hne), 
X2 ^ Xil^l~ (dashed hne), Li^t = 10^ pb^^ 
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Figure 6.5: Domains where the observed edge in the Mi+i-^ distribution is due to the 
decays X2 ~^ ^l^"^ ~^ Xil~^l~ (sohd hue), ~^ ^r^^ ~^ Xi^^^~ (dashed-dotted hue), 
X2 ~^ Xi^^^~ (dashed hne), L^t = 10"^ pb~^. 



74 




50 100 150 200 250 300 350 400 450 500 



mo (GeV) 

Figure 6.6: Domains where the observed edge in the M;+;- distribution is due to the 
decays X2 ~^ ^l^"^ ~^ Xil~^l~ (sohd hue), ~^ ^r^^ ~^ Xi^^^~ (dashed-dotted hue), 
X2 ~^ xV'^l' (dashed hue), Li^t = 10^ pb"^ 
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Figure 6.7: Expected Standard Model background in inclusive 3 lelpton final state. 
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Figure 6.8: Expected mass spectrum for mSUGRA Point 4: mo = 200 GeV, mi/2 = 
100 GcV, tan/9=2, = and /j, < 0. a) Leptons with pr > 15 GeV in \r]\ < 2.4 are 
considered; b) as in previous case, but two leptons which enter invariant mass spectrum 
are required to be isolated. The dotted histogram corresponds to the SM background. 
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Figure 6.9: Expected mass spectrum for mSUGRA Point 1: mo = 100 GeV, 
mi/2 = 300 GeV, tan/?=2.1, = 300 GeV and > 0. The dotted histogram corresponds 
to the SM background. 
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Figure 6.10: Expected /"''/" mass spectrum for mSUGRA Point 5: niQ = 800 GcV, 
mi/2 = 200 GeV, tan/3=10, Aq = and > 0. The dotted histogram corresponds to the 
SM background. 
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Figure 6.11: Dilepton mass spectrum for mSUGRA Point from different (w^ioj ^^^1/2) 
regions in inclusive 3 leptons and 3 leptons + -E™"'' cliannels. Contribution from SM back- 
ground is shown witli daslied liistogram. Ttie otlier mSUGRA parameters are: tan/3=2, 
Ao = and < 0. 
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mSUGRA parameters: tan^=2, Ar.-0, )Ll<0 




Figure 6.12: Explorable region of mSUGRA in inclusive 3 leptons and 3 leptons + 
j^miss Q];ia]2nels at Lint = 10"^ pb~^ and Lint = 10^ pb~^, respectively. Shaded regions are 
excluded by theory (mo ^ 80 GeV) and experiment (mi/2 ^ 90 GeV). The cosmologically 
preferable region 0.15 < Qh"^ < 0.4 is also given. Simulated mSUGRA Points are shown 
as well. 
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7 Exclusive 2 leptons + no jets + Ej^'^^^ channel — 
search for sleptons 

The previous two sections described inclusive searches suitable for the discovery of SUSY. 
We now focus on more exclusive final states where individual SUSY particles may be 
targeted. 

The slepton production and decay features described in section 2.2 give rise to several 
interesting experimental signatures. In this study we concentrate on the possibilities 
to detect a slepton signal in the simplest event topology two leptons + no jets 

taking into account all possible SUSY and SM contributions to this topology 0. The 
observability of the slepton signal is evaluated summing all flavors and all production 
mechanisms which contribute to the two same-flavor opposite- sign leptons + EJ^^'^'^ + no 
jets final state. Evidence for a slepton signal will be judged from an excess over full 
SUSY and SM backgrounds, but presence of a SU SY signal will be estimated from an 
excess over only the expected SM backgrounds. 

The main expected Standard Model background is tt production, with both W^s de- 
caying to leptons, or one of the leptons coming from a W decay and the other from the 
6-decay of the same t-quark. The other SM backgrounds considered are Wtb, WW, WZ, bb 
and rr-pair production, with decays to electrons and muons. The assumed cross-sections 
are: att=660 pb (M^ = 175 GeV), crvi/fb=160 pb, aww=7^ pb, o"vfz=26 pb, (T^t-=7.5 pb, 
(7^5=5000 pb (p^ > 100 GeV). No forcing of decay modes is implemented. 

As internal SUSY backgrounds we consider processes qq, gq, gg which, through cascade 
decays with jets outside the acceptance, below the detection threshold, or lost in cracks, 
can also lead to two leptons -\- no jets final states. The xfx2 process for example 

leads mainly to three-lepton final states. However, it can also contribute to two-lepton 
final states if one of the leptons is outside the acceptance or is misidentified, or if only 
the X2 decays to leptons, and jets from the xf decay are not reconstructed. 

We perform a mapping of the mSUGRA parameter space to determine the optimal 
set of signal selection cuts in the various regions of (mo,mi/2) plane, taking into account 
event kinematical properties, and evaluate the observability of the signal. An important 
issue is to find out if there is a gap in slepton detectability, especially in the low mass range 
at the junction between the observability domains at LEP2, Tevatron and the LHC. The 
mSUGRA-parameter space points (mo,mi/2) probed are shown in Fig. 7.1. This figure 
also shows the masses mj^ and m^o at each point, the theoretically and experimentally 
excluded regions and the explorable domains in sparticle searches at LEP2 (300 pb^^) 
and the Tevatron (1 fb~^). 

Taking into account the signal and background event topologies the following kine- 
matical variables are found useful for slepton signal extraction: 

i) for leptons: 

•pip-cnt on leptons and lepton isolation {I sol), which is defined here as the calori- 
metric energy flow around the lepton in a cone of AR < 0.5 divided by the lepton energy; 

• effective mass of two same-flavor opposite-sign leptons, to suppress WZ and 
potential ZZ backgrounds by rejecting events in a mz ± SMz band; 

• A0/+/-, the relative azimuthal angle between two same-flavor opposite-sign lep- 
tons; 
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ii) for Ef^'': 

•A</)(£^™*^, /"•"/"), the relative azimuthal angle between EJ^'^'^'^ and vector sum of 
the transverse momenta of the two leptons. 

hi) for jets: 

• "jet veto"-cut: no jet with Ej^^ lower than certain threshold, in some rapidity 
interval, typically | rij^t \< 4.5. 

The data selection cuts based on these variables are separately optimised in various 
parts of the explored parameter space For illustration we show some typical distribu- 
tions. Figure 7.2 shows the E'™'^^ distributions for signal and various SUSY backgrounds 
at two (mo, rni/2) points, one with low and one with high slepton masses. The hardest 
spectra in both figures correspond to qq production and the softest ones to chargino- 
neutralino production. This figure also shows the main SM backgrounds. For low slepton 
mass there is no big difference in shape between signal and SM-background spectra, thus 
it is not very advantageous to apply a hard i?™***^ cut. 

Figure 7.3 shows the A(f){EJp^^^ distribution at the point from the low slepton 
mass region, ttiq = 86 GeV, mi/2 = 85 GeV, after p^'^ , cuts on lepton isolation and E^""^^. 
The internal SUSY and SM backgrounds are still much larger in magnitude than the 
signal. Keeping only events with large relative azimuthal angles significantly improves 
the signal to background ratio. 

Figure 7.4 illustrates the expected rejection factors which can be obtained by a "cen- 
tral" jet veto for two E^^^ thresholds as a function of jet acceptance at point mg = 86 GeV, 
mi/2 = 85 GeV. Signal acceptance decreases slightly with increasing \rijet\ of the jet veto 
region, whilst for sufficient rejection of SUSY and SM backgrounds {tt) it is important to 
have calorimetric coverage up to \ rj \= 3.5 — 4.0. The rejection power against g, q and ti 
improves by a factor of 1.5-1.8 with increasing coverage from | |< 2.5 to | |< 3.5. To 
achieve significant background rejection reliable jet detection down to E^ ~ 30 GeV is 
also essential. 

Figure 7.5 shows the Acpi+i- distribution for the main SM and SUSY backgrounds 
and events with direct slepton production at the same point mo = 86 GeV, mi/2 = 85 
GeV. For the signal, small A0;+;- events come from uiui production whilst events with 
large A0;+;- are due to charged slepton production. The main contribution to the slepton 
signal from X2 production in the two- lepton final state comes from decays of thus 
the distribution of these events also peaks at low relative azimuthal angles. All processes 
which go through X2 decays to leptons, including g, q contributions have almost the same 
shape A(^-distribution. In some parts of parameter space it may be more advantageous to 
keep only events with large A0;+;-, which rejects most of the internal SUSY backgrounds. 
However, this cut also eliminates direct z> production events and the indirect contribution 
from X2 decays. Only xfxf events with leptons coming from different parent charginos 
survive as a significant SUSY background. 

To extract the slepton signal at points from domain I in Fig. 2.8b (mo ^ 0.5mi/2) 
with mj^ ~ 100 GeV, just beyond of the LEP2 and Tevatron sparticle reaches, where the 
SM backgrounds are dominant, the selection Set is appropriate at low luminosity. The 
Set 0' is better adapted for the points from domain II (mo ^ 0.5mi/2), see Table 7.1. 

The expected number of signal and background events in this particular difficult region, 
after all cuts, is shown in Table 7.2. We also give the signal to background ratio and 
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Table 7.1. Low mass slepton selection criteria 
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^<pm- > n 


130 


no cut 



significance for slepton production and the same quantities for overall SUSY production. 

It is clear that in this domain of parameter space, the slepton (m^- ~ 100 GeV) con- 
tribution to the finale state topology is only a modest one. The slepton signal can be 
detected from an excess of events over background expectations, but not from any par- 
ticular kinematical feature. The S/B is about 0.3. Of the SM backgrounds, WW is the 
largest and essentially irreducible, tt and MVtb are comparable and at the same level as the 
signal after hard jet veto cuts. The main SUSY background is due to xf • the point 
mo =85 GeV, mi/2 =86 GeV, from domain I, the total number of signal events iV^'"*"* is 
just due to direct slepton production, i.e. there are no indirect slepton contributions, but 
there are additional SUSY contributions to this topology. 

For points with mo ^ 0.5mi/2 (domain II), no A0;+;- cut is applied so as to take 
advantage of all SUSY processes contributing to slepton production. iV|'^~*°* in this case 
has contributions from both direct and indirect slepton production. The contribution 
from indirect slepton production is comparable to the direct one and helps to reach a 5a 
significance level needed for slepton signal visibility. In this region the only backgrounds 
are SM contributions. The irreducible WW is still the dominant, and Wtb increases 
significantly. 

To achieve the highest mass reach with 10'^ pb^^ a somewhat harder set of cuts is 
appropriate, (Set 1 in Tab. 7.3). The reach in parameter space with W^pb^^ is shown 
in Fig. 7.6; it corresponds to Icft-slcpton masses up to ~ 160 — 180 GeV, i.e. well above 
the LEP2 reach {mf ^ 95 GeV). In the region mo ^ 0.5mi/2 (domain II) the total slepton 
signal coincides with the overall SUSY signal, but for points from region mo ^ 0.5mi/2 
(domain I) the additional SUSY processes contribute to the overall SUSY signal, but not 
to the slepton signal. 

With 10^ pb^^, a significantly higher slepton mass range becomes accessible. The 
average lepton transverse momentum and missing energy increasing with m^- become sig- 
nificantly harder than for the main backgrounds. Harder kinematical cuts are therefore 
preferable. Slepton production becomes the main and ultimately the dominant contribu- 
tion to this topology. We use four main different sets of cuts (Table 7.3) optimised for 
different regions of (mo,mi/2) parameters space. These sets are chosen taking into ac- 
count the mass relations, production mechanisms and decay patterns discussed in section 
2.2. 

"Set 2" is used to explore the mf ~ 200 GeV domain. The ti, Wtb backgrounds are 
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Table 7.2: Expected number of direct (A^f"'^*'^) and indirect (ATf'"^'^^*'-) slepton signal 
events, and SM (A^|*0 and SUSY background (N^^'^^) events, after cuts, at 10^ pb'^ 

The signal to background ratio A^|^'~*°*/^b* and significance of slepton signal 
^si-tot I ^ ^si-tot ^ j^tot ^igo gjygj^. ^i^g g^^g for overall SUSY signal Ar|^^^/Ar|Af^ 



(mo,mi/2) ^ 


(86,85) 


(4,146) 


(20,160) 


(20,190) 




Set 


Set 0' 


Set 0' 


Set 0' 


// 


323 


315 


319 


214 




14 


45 


48 


40 




46 


253 


217 


168 




13 


8 


— 


— 




35 


6 


— 


— 







2 


— 


— 


WW 


454 


1212 


1212 


1212 


Wt6 


163 


577 


577 


577 


tt 


345 


574 


574 


574 


wz 


15 


A O 

43 


43 


A O 

43 


TT 


15 


15 


15 


15 




323 
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214 
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314 


265 


208 




431 
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108 










992 


2421 
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0.3 


0.26 


0.24 


0.17 


^Sl-tot 1 ^ ^Sl-tot ^ ^tot 


8.6 


11.4 


10.7 


7.9 




0.4 


0.26 


0.24 


0.17 




11.4 


11.4 


10.7 


7.9 



suppressed with hard cuts on pt of the leptons, i?™*^ and the "jet veto". This last cut 
eliminates entirely the backgrounds from strongly interacting sparticles. We keep events 
with small relative azimuthal angle between the leptons mainly to improve the 

rejection of SM backgrounds, and to eliminate xtxi production when it is necessary. The 
slepton S/B ratio is about 1. 

"Set 3" is used for m^-^ ~ 300 GeV at points with mo ^ 0.5mi/2, whilst for points 
with rriQ < 0.5mi/2 "Set 4" is more appropriate, as the cross-sections are smaller. Fur- 
thermore, there is no source of indirect slepton production. "Set 4" is also used to 
explore the slepton mass range 300 GeV ^ m^- < 350 GeV. "Set 4"' is used for masses 
350 GeV < mj^ < 440 GeV, and "Set 4"" for^m,-^ > 440 GeV. Details of the signal and 
background contribution at each simulated point can be found in 0. 

Figure 7.6 summarizes the results of this study. At each point investigated the expected 
slepton signal significance, na, and the signal to background ratio (in parentheses) is given. 
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Table 7.3. High mass slepton selection criteria 





Set 1 


Set 2 


Set 3 


Set 4 


Set 4' 


Set 4" 


Isol < 


0.1 


0.1 


0.1 


0.1 


0.1 


0.03 


p^'^ > (GeV) 


20 


50 


50 


60 


60 


60 


Ejntss^ (GeV) 


100 


100 


120 


150 


150 


150 




150 


150 


150 


150 


150 


150 


No jets in \r]jet\ < 4.5 














with E^"* > (GeV) 


30 


30 


30 


45 


45 


45 


No M;^ eMz±5m (GeV) 


91 ±5 


91 ±5 


91 ± 5 


91 ±5 


91 ±5 


91 ±5 
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130 


130 
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for an integrated luminosity of 10 pb^ . The 5a reach is shown as well. Recall that the 
selection procedure has been optimised in terms of signal significance Ns / a/ Ns + Nb to 
investigate the maximum slepton mass reach, and not in terms of sample purity (i.e. 
S/B ratio). The reason is that if SUSY were found at Fermilab where g/q masses up to 
< 400 — 450 GeV can be probed, sleptons can be probed at the Tevatron and LEP2 only 
up to a mass of mj < 100 GeV, thus the main task of LHC would be to explore slepton 
spectroscopy, with masses up to riig/g < 400 GeV. Figure 7.7 shows the same information 
as Fig. 7.6, but in terms of left slepton and lightest neutralino LSP masses. The domain 
that can be explored extends up to m^-^ ~ 340 GeV for all allowed misp {< 200 GeV), 
and up to mj^ ~ 340 — 440 GeV only if ttilsp ~ (0.45 — 0.6) ■ rrij^ for a given ttlj^. The 
mass reach for right sleptons is up to 340 GeV. As the main SM background is WW 
production, to study the stability of our results we assumed a 50% uncertainty in the 
WW cross-section. Figure 7.8 shows how much the 5a reach boundary would be modified 
if the WW cross-section were twice as large as predicted by PYTHIA. 

A very interesting aspect of slepton searches at LHC from the cosmological point of 
view is that, in this topologically simplest channel, the explorable domain for 10^ pb~^ 
almost coincides with the regions of parameter space where the LSP would be cosmo- 
logically relevant (Figs. 7.6 and 7.7). The boundaries for the domains we can probe in 
slepton searches and the neutralino dark matter Qh"^ density contours have shapes similar 
to slepton isomass curves as both are largely determined by slepton masses. In particular, 
neutralino annihilation is governed by slepton exchange with aann ~ 1/mj. Too large a 
slepton mass would induce too small an annihilation cross-section and thus too large a 
relic density which would overdose the Universe. The Qh"^ range prefered by mixed dark 
matter scenarios is for 0.15 < Qh'^ < 0.4 (with the Hubble scaling constant 0.5 ^ h? ^ 



0.8) which corresponds to a slepton mass range 100 GeV ^ mj < 250 GeV, (see |]TT|, [p7| .) 

We conclude that to search for direct production of sleptons the most appropriate 
channel is 2 leptons+ EJj^'^'^^ and no jets. In regions of low mo,mi/2, just beyond the 
sparticle reaches at LEP2 and the Tevatron, the slepton signal can be detected, but the 
search will be difficult as I is not the dominant contribution to this final state. In signifi- 
cant parts of the explorable (mo, ^1/2) space, however, direct slepton pair production is 
the main contribution to this final state, with non- negligible contributions from indirect 
production too, mostly through xfjX^- The inclusion of this indirect slepton production 
improves the mass reach. With 10^p6~^ luminosity, CMS is sensitive up to mi < 160 
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GeV. With 10^ pb^^ luminosity the reach extends up to m^-^ ~ 340 GeV for all allowed 
rriLSP (< 200 GeV), and up to m^-^ ~ 340-440 GeV if rriLSP ~ (0.45-0.6) -m;-^ for a given 
m^-^. Throughout the domain where the slepton signal should be detectable, it is always 
contaminated by backgrounds. The dominant background is always WW production. In 
the mf^ ~ 100 GeV domain the WW background exceeds the number of signal events, 
whilst for mj^ > 200 GeV the number of background events is comparable or smaller than 
the expected number of signal events. The "jet veto" cut effectively eliminates the g, q 
backgrounds, but also the potentially important indirect slepton production from g and 
q. This issue requires a separate study. 

In the region mo ^ 0.45 •mi/2 sleptons decay directly into the LSP and sneutrinos thus 
decay invisibly, z/ — > Xi^- this portion of parameter space only charged sleptons can 
be detected. As right sleptons are also significantly lighter than left ones they are mostly 
eliminated by the hard selection cuts needed to suppress the backgrounds, so mainly left 
sleptons are present in the final event sample. Left sleptons will thus be detectable in the 
entire region where a slepton signal is seen, whilst right sleptons will be detectable mainly 
in domain nio > 0.45 • mi/2. Sneutrinos will contribute to this final state only in a very 
limited area with mi/2 ^160 GeV and nio < 150 GeV. 

The search for sleptons is important and significant, as LEP2 and the Tevatron can 
only explore the slepton mass range up to m^-^ ~ 100 GeV. The domain of (mo, mi/2) 
parameter space explorable through this simplest experimental channel, with two same- 
flavor, opposite-sign leptons -h EJp"^^^ + no jets, covers almost the entire domain of param- 
eter space where SUSY is plausible in general terms (minimal fine-tuning) and covers all 
of the parameter space where SUSY would be of cosmological significance (at least for 
tanl3 = 2, Ao = 0, /X < 0), with > lO^pft"-*^ of integrated luminosity. 
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Slepton mapping of parameter space 



mSUGRA-MSSM ; tg p= 2, Aq = 0, 



600 



500 



400 



> 

^ 300 

e" 

200 



100 



/ / / 1 \ 455,239 
/ / A \ • 


/ / /] 406,216 •V41 6,21 6 — 




M *+-Jt-, I C7H- 


A • 314,181 \ 
/ / / 378,1 72 • 


• 408,172 _ 




Qh^= 1 

K 


7^^/^212,123 

/ / #212,107 i 

// |_^p2Tevatron / 

-7143,82 / 7^212,44 / 

/*1 23,69 / 1/ / 
/ 0212,69/ 

64 ^ - _ ^, V " - 


<► 344,135 _ 
(m^l^,m~o) at 
points investigated 

• 314,91 — 

311,44 344,44 
• • 















100 



200 300 
mo (GeV) 



400 500 



Figure 7.1: Points in mSUGRA parameter space where slepton production has been 
studied. The numbers at each point give corresponding values of m^-^ and m^o^. The 
shaded areas represent excluded regions on the bases of LEP and Tevatron data. The 
dashed lines represent the expected reach with sparticle searches of LEP2 and Tevatron 
with 1 fb~^. The VLh? = 1 absolute upper limit contour for neutralino dark matter is also 
shown. 
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distributions at points: 
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Figure 7.2: i?™*'^ distribution for slepton signal (DY production), SUSY backgrounds 
{9i 'X.f'X.2) parameter space points: a) uiq = 86 GeV, mi/2 = 85 GeV, b) uiq = 97 

GeV, mi/2 = 420 GeV and c) the main SM backgrounds {WW, tt). 
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Figure 7.3: Relative azimuthal angle between E^^^^ and the resulting dilepton mo- 
mentum for the slepton signal and main SM and SUSY backgrounds in 2leptons + E^^^^ 
+ nojets final states. 
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Figure 7.4: Rejection factors a) for internal SUSY backgrounds, b) for SM backgrounds 
by central "jet veto", and c) slepton signal acceptance, as a function of jet rapidity 
acceptance | rjjet \ and for two jet detection thresholds. 
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Figure 7.5: Relative azimuthal angle between two same-flavor opposite-sign leptons at 
one point on parameter space, a) for the various contributions to the SUSY slepton signal 
and for a SUSY background {xfxt)^ b) for the main SM backgrounds, in 2leptons + 
j^mtss _|_ jg-^^ ^nal states. 
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Slepton mapping of parameter space 
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Figure 7.6: Expected reach in parameter space for slepton searches in the 2leptons+ 
E^^^'^+no jets final states with lO'^pfe"^ and 10^ pb~^ in CMS. The reach is given as a 5cr 
contour. At each point investigated the slepton signal significance is indicated as na, and 
the S/B ratio is given in parenthesis, where B is the full SM + SUSY background to the 
slepton signal. The cosmological absolute upper bound (fi^o = l, h=l) density contour for 
neutralino dark matter Qh"^ = 1 is indicated. 
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Explorable range masses 
in slepton searches in CMS 
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Figure 7.7: Same mapping as in Fig. 7.6, but in terms of masses of the left slepton 
and the LSP. 
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Slepton mapping of parameter space 
mSUGRA-MSSM, tan p= 2 Ao= 0; |i< 
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Figure 7.8: Modification of tlie 5a reacli in slepton searclies if the production cross- 
section for the main SM background of WW pair production is increased by a factor of 
2. 
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8 Exclusive 3 leptons + no jets + Ej^'^^^ channel — 
search for chargino/neutralino pair production 

Direct production of xf X2 followed by leptonic decays of both sparticles (see section 2.1) 
gives three high-p^- isolated leptons accompanied by missing energy. These events are free 
from jet activity except for those coming from initial-state QCD radiation. Therefore we 
search for X1X2 P^^^ production in 3 leptons + no jets+ i?™'^* final states p. 

We have considered the following SM processes contributing to the background: WZ, 
ZZ, ti, Wtb, Zbb, bb with the cross-sections awz = 26 pb, azz = 15 pb, au = 670 pb for 
a top mass of 170 GeV, awtb = 160 pb, azbi = 600 pb, cr^^ = 500 /ib. Three leptons in 
the final state are obtained by folding the appropriate leptonic decay modes. We have 
forced Z W^ l^v and b — civ decays, where / = /x, e, r. In case of ti and 

Wife only one W has been forced. These SM backgrounds, except ZZ production, can be 
suppressed by requiring that all three leptons be isolated, and vetoing central jets. ZZ 
events contain four isolated leptons and asking for exactly three isolated leptons in the 
event leaves only a small fraction of the ZZ continuum, when one of the leptons escapes 
detector acceptance. The WZ, ZZ and Zbb backgrounds can be suppressed by removing 
events with a dilepton invariant mass close to the Z mass. At each investigated point 
in the (mo,mi/2) plane we have studied all SUSY production processes which can lead 
ultimately to three leptons in the final state. These are the strong interaction production 
of gg, gq, qq pairs; associated production of gx, qx pairs; chargino-neutralino XX pair 
production, and slepton pair production, //, W, w. The SUSY event sample with three 
high-pr isolated leptons is dominated by strong production; but, of course, the jet veto 
requirement is very efficient in reducing g and q events which in their cascade decays 
produce many central high-i?^ jets. It has to be pointed out that SUSY events other than 
the directly produced xf X2^iid which survive all selection criteria, may still contain some 
xf ^2'^ produced in cascade decays of heavier sparticles. These events can be considered 
as an "indirect" xf X2 signal, a situation similar to the one encountered in discussing 
slepton searches. 

In Fig. 8.1 we show the fully simulated SUSY points in the (mo,mi/2) plane, as 
representatives of various regions according to the predominant or characteristic decay 
modes. In order to illustrate event kinematics we choose three points from different 
(mo,mi/2) regions among those shown in Fig. 8.1: 

51 (mo, mi/2) = (100 GeV, 100 GeV) with 
M(xf ) ^ M(x^) = 98 GeV, M(x?) = 45 GeV, 
M(^) = 294 GeV, M(g) = 270 GeV 

52 (mo,mi/2) = (2000 GeV, 175 GeV) with 
M(xf ) ^ M{xl) = 154 GeV, M(x?) = 76 GeV, 
M(^) = 551 GeV, M(g) = 2032 GeV 

53 (mo,mi/2) = (100 GeV, 400 GeV) with 

M(xf ) ^ M(x^) = 341 GeV, M(x?) = 172 GeV, 

M(^) = 983 GeV, M(g) = 861 GeV 
Figure 8.2 shows the da/ dpip distribution in the three-lepton event sample with pip > 15 
GeV. Distributions are shown for signal and SM background. At this stage the dominant 
SM backgrounds are ti and Zbb. For high mi/2, e.g. point S3, the lepton pt cut can be 
set at 20 GeV. The loss in signal acceptance is ~12% compared to a 15 GeV cut, whilst 
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tt and Zbb are reduced by an additional factor of ~2. 

To suppress background events with leptons originating from semileptonic decays 
of 6-quarks, we apply tracker isolation, requiring all three leptons to be isolated. If 
there is no charged track with greater than a certain threshold p^* in a cone of 
R = Ar]'^ + = 0.3 about the lepton direction, it is considered as isolated. Figure 
8.3 a shows the lepton isolation rejection factor against ti background versus the 
on charged tracks. The two curves correspond to two different instantaneous luminosity 
regimes, L = 10^^ cm~^s~^ and L = 10^^ cm~^s~^. In the latter case we have superimposed 
15 "hard" pile-up events. Figure 8.3 b shows signal acceptance versus p^*. We choose p^* 
= 1.5 GeV and 2 GeV for low and high luminosity running, respectively. The correspond- 
ing rejection factors against ti are about 60 and 50, whereas signal acceptance decreases 
from about 90% to 70%. The next requirement is the central jet veto, which allows us 
reject the internal SUSY background coming from g and q cascade decays overwhelming 
xfx2 direct production. This cut also reduces some SM backgrounds. In Fig. 8.4 we show 
jet multiplicity and jet Et distributions for signal, SUSY gg, gq, qq background events at 
point SI, and for the ti background. Since at this {1^0,7711/2) point the gluino is heavier 
than the squark, it decays to a squark-quark pair. Therefore gluino events are richer in 
jets compared to squark events. Figure 8.5 shows jet veto rejection factors against gg, 
gq, qq and ti backgrounds and the signal acceptance as a function of jet pseudorapidity, 
|i^jet|^ for two jet detection thresholds: Ej^* > 25 GeV and i?jf* > 30 GeV. As the jet veto 
region is decreased down to \ri^'^^\ ^ 3.5, the rejections/acceptance are not too sensitive 
to jet rapidity coverage, but below this value of \r]^'^^\, rejection factors decrease rapidly, 
degrading significantly the separation of direct X1X2 production from internal SUSY and 
SM backgrounds. The reconstruction of jets at \r]^^*\ ~ 3.5 implies calorimetric coverage 
up to |?7| ~ 4.3 at least, as the jets are collected in a cone of AR = 0.8. With increasing 
mo and/or mi/2 squarks and/or gluinos become heavier. Thus jets are harder and SUSY 
events are easier to suppress by a jet veto with relaxed requirement on \ri^^^\ and E^^^. 
Furthermore, for the signal events the QCD radiation becomes harder {xf, X2 more 
massive) and the low production cross-section requires the highest LHC luminosity, with 
an event pile-up contaminated environment. For this region, the jet veto threshold can be 
set with Eq^^ as high as 50 GeV and is justified by the need to suppress the SM (mainly 
ti) background. 

Despite the fact that missing energy is always present in SUSY signals (at least in 
i?-parity conserving models), for some regions of mSUGRA parameter space, it is advan- 
tageous not to apply any specific cut on it. Figure 8.6 shows the distribution of ii^™*^ in 
events with three isolated high-p^^ leptons and no jets for the signal and the two main SM 
backgrounds: ti and WZ. For point SI the shape of the distribution is very similar to 
that of the background and therefore no cut on El^^^^ is applied. For point S3 requiring 
j^miss-^ 80 GeV significantly reduces ti and W Z backgrounds whilst keeping more than 
75% of signal events [§. 

For WZ background suppression, which after all the previous requirements remains 
the dominant SM background, an Mz cut is necessary. In the signal events, the shape of 
the (same-flavor opposite-sign) dilepton invariant mass Mi+i- distribution is determined 
by the X2 decay kinematics and has a very sharp edge. Figure 8.7 shows the Mi+i- 
distribution for the signal at three parameter space points SI, S2, S3 and for WZ events. 
In events having all leptons of the same flavor, the combinatorial background obscures 
the sharp mass edge in the signal events. A rejection factor of ~ 17 for WZ is obtained by 
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Table 8.1 


: Selection criteria 


Set 1 


Pt cut 


3/ with pir > 15 GeV in |V| < 2.4(2.5) 


appropriate for 


Isolation 


No track with pt > 1.5 GeV in R = 0.3 


mi/2 < 140 GeV 


Jet veto 


No jet with E^^' > 25 GeV in < 3.5 




Z-mass cut 


Mi+i- < 81 GeV 


Set 2 


Pt cut 


3/ with > 15 GeV in |V| < 2.4(2.5) 


appropriate for 


Isolation 


No track with pt > 1.5 GeV in R = 0.3 


140 GeV < mi/2 


Jet veto 


No jet with ' > 30 GeV in |r/^'^*| < 3 


< 180 GeV 


Z-mass cut 


Mi+i- < 81 GeV 


Set 3 


Pt cut 


3/ with pi^ > 15 GeV in |V| < 2.4(2.5) 


appropriate for 


Isolation 


No track with pt > 2 GeV in R = 0.3 


180 GeV < mi/2 


Jet veto 


No jet with E^^' > 40 GeV in \ri^^^\ < 3 


< 300 GeV 


missing Et 


E^'''> 50 GeV 




Z-mass cut 


Mz ± 10 GeV 


Set 4 


Pt cut 


3/ with pi^ > 20 GeV in |r/'| < 2.4(2.5) 


appropriate for 


Isolation 


No track with Pt > 2 GeV in R = 0.3 


mi/2 > 300 GeV 


Jet veto 


No jet with E^^* > 50 GeV in < 3 




missing Et 


EJp'''> 80 GeV 




Z-mass cut 


Mz ± 10 GeV 



removing events with M;+/- > Mz — 10 GeV, whilst the small loss of signal events, < 10% 
for point SI, comes from events when all three leptons have the same flavor. However, 
for higher values of mi/2, e.g. point S3, the endpoint of the Mi+i- distribution is above 
the Z peak; in this case we cut only events falling within the Z mass window: Mz — 10 
GeV < Mi+i- <Mz + 10 GeV. 

On the basis of the above considerations we adopt four sets of selection criteria, ap- 
propriate for different (mo,mi/2) regions, given in Table 8.1. 

The mi/2 boundaries are approximate and, for optimal cuts, they depend also on mo. 

Except for Set 4, which we apply for rather massive xfx2^ the pt cuts on leptons are 
the same, i.e. pip > 15 GeV. Set 1 and Set 2 are appropriate for the regions where the 
signal is accessible already at low luminosity. In the regions where we apply the cuts of 
Set 3 or Set 4 a high luminosity is required, as they correspond to smallest cross-sections; 
the luminosity determines the isolation requirement. With increasing mi/2 the jet veto is 
relaxed, varying the threshold from 25 to 50 GeV. A cut on EJp'^^^ is applied only for 
Set 3 and Set 4; and finally, the Mz cut is looser for higher mi/2. 

Table 8.2 gives the SM background cross-sections for these selection criteria after each 
cut. The two main SM backgrounds surviving these cuts are WZ and tt. With Set 1 and 
Set 2 the dominant background is WZ and with Set 3 and Set 4 it is tt. 

Details on each simulated point can be found in ||^ . The signal significance is defined 
as the number of signal events divided by the square root of all events (SUSY + SM) 
passing our selection criteria. With this definition we give: i) evidence for existence of 
Supersymmetry in 3/ + no jets+ E^^^^ channel; ii) overall evidence for X1X2 production 
and iii) evidence for direct production of xfx2- The ratio of direct xfx2 events over other 
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Table 8.2: SM background cross-sections (in fb) at each step of applied cuts. 



Set of cuts 


Uuts 


W Zj 




tt 


WW 


ZiOO 


00 


Set 1 


Vt 


111 


21.7 


2544 


139.1 


731 


9300 


Isolation 


102 


19.7 


39.8 


4.09 


28.8 


< 0.30 


Jet veto 


73.8 


11.6 


4.12 


1.71 


14.6 


< 0.15 


Mz 


4.iy 


U.oy 


QIC 


U.oD 


U. / D 


< U.io 


Set 2 


1 

Ft 


111 


21.7 


2544 


139.1 


731 


9300 


Isolation 


102 


19.7 


39.8 


4.09 


28.8 


< 0.30 


Jet veto 


76.1 


12.7 


4.84 


1.90 


16.1 


< 0.15 


71 T 

Mz 


4.49 


0.99 


2.48 


1.05 


0.91 


< 0.15 


Set 3 


Pt 


111 


21.7 


2544 


139.1 


731 


9300 


Isolation 


75.4 


14.7 


46.3 


4.82 


32.9 


< 0.45 


Jet veto 


51.3 


8.69 


4.77 


1.91 


15.7 


< 0.40 


rpTTits s 

Ulrp 


22.0 


2.61 


3.01 


0.48 


1.42 


< 0.01 


Mz 


1.93 


0.35 


2.63 


0.40 


0.06 


< 0.01 


Set 4 


Pt 


88.4 


16.6 


1330 


65.2 


378 


2800 


Isolation 


60.5 


11.4 


19.0 


1.33 


10.3 


< 0.04 


Jet \-eto 


16.1 


7.92 


3.25 


0.53 


1.78 


< 0.01 


rpmiss 

Ulrp 


5.77 


0.68 


0.73 


0.05 


0.30 


< 0.001 


Mz 


0.45 


0.09 


0.60 


< 0.05 


< 0.04 


< 0.001 



events surviving the cuts, N-±.q/[Nall- N-±-o), is also given. All investigated points are 

Xi X2 Xi X2 

analyzed in a similar manner and 5a significance contours are obtained in the (mo,mi/2) 
plane, as shown in Fig. 8.8. 

In the regions of small mo and mi/2 200 GeV and < 150 GeV, respectively) indi- 
rectly produced xfx2 events contribute significantly in the 3/ + no jets+ EJp'^^^ channel. 
However, this contribution is expected to be less than ~15% of the overall xfx2 Produc- 
tion. Therefore, we can give 5a significance contours for direct X1X2 production. 

Prom Fig. 8.8 one can see that at a low integrated luminosity of L^nt — 10^ pb~^ the 
direct production of xfx2 can be distinguished from background up to mi/2 ^150 GeV 
for all rriQ. A further increase of luminosity up to Lint = 10^ pb~^ extends the explorable 
region only by about 10 ^ 20 GeV for mo ^ 120 GeV. However, for mo ^ 120 GeV the 
gain in parameter space is much larger - up to mi/2 ^ 420 GeV, reflecting the more 
advantageous behavior of the production cross-section times branching ratio shown in 
Fig. 2.4. 

It should be pointed out that the internal SUSY backgrounds are particularly dan- 
gerous for small mo, mi/2, in the region mo ^ 200 GeV, mi/2 ^150 GeV, generating a 
background higher or comparable to the SM contribution because of the high production 
cross-sections of the relatively light gluinos and squarks. Softness of these gluino/squark 
events makes it difficult to suppress them below the SM level, even with a very strict jet 
veto requirement (selection criteria Set 1). For some regions of the (mo,mi/2) parameter 
plane, chargino/neutralino and slepton pair productions are comparable to gluino/squark 
contributions. As mo and mi/2 increase, more massive gluinos and squarks result in 
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harder jets, and decreased cross- sect ions, and hence the SUSY background becomes neg- 
hgible. At the 5a boundary with an integrated luminosity of 10^ pb~^, background events 
dominantly consist of SM processes. At this boundary the signal to background ratio 
is expected to be > 0.20 with about 150 signal events {Lint = 10^ pb~^) for mo ^ 120 
GeV. In this region the dominant SM background is WZ production which cannot be 
suppressed by a stronger jet veto or lepton isolation requirements. For mo ^ 120 GeV 
the signal to background ratio is > 0.60 with about 70 signal events. Here the dominant 
SM background is ti production. The lower number of signal events needed to attain a 
5a significance and the higher signal to background ratio in this region is the result of the 
harder event kinematics of xtx2^ allowing for more stringent selection criteria (Set 4). As 
a consequence the SM contamination is smaller 0. 

At small mo, the xf X2 reach in parameter space is very luminosity-dependent. This 
region of small mo is of particular importance, since it is the prefered one from the 
cosmological point of view, where x? is the candidate for the cold dark matter in the 
Universe. With the present lifetime of the Universe (> 10^° years) and the estimate 
of the relic neutralino abundance at freeze-out time as a function of (mo, Tni/2), one 
can determine regions of parameter space where Qh'^ < 1; here Q is the ratio of the 
relic neutralino density to the critical density which closes the Universe and h is the 
Hubble constant scaling factor, 0.5 < h < 0.8. In Fig. 8.8, the fi/i^ = 0.15 and flh'^ = 
0.4 contours mark the most probable region of the parameter space, corresponding to a 
mixed dark matter scenario with 60% of neutralino dark matter. One can see that xf'X.2 
direct production covers a significant part of this prefered region of parameter space. 

We can conclude, that: 

• by measuring the excess of events in the 3/ + no jets+ E'^*'''^ channel over SM and 
internal SUSY backgrounds, X1X2 direct production can be investigated in mSUGRA 
parameter space up to mi/2 ^ 170 GeV (150 GeV) for all mo at a luminosity L^nt = 10^ 
pb-i (10^ pb-^); 

• with a high integrated luminosity of Lint = 10^ pb^^ and for mo ^ 120 GeV, the 
discovery region for direct production of xf X2 extends up to mi/2 ^ 420 GeV. 
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Figure 8.1: Simulated SUSY points for chargino/neutralino pair production. Sigma 
times branching ratio contours are also shown. Shaded regions are excluded by theory 
(mo ^ 80 GeV) and experiment (mi/2 ^ 90 GeV). 
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Figure 8.2: dajd-p^j, distributions for signal and SM background. Left plots show- 
distributions for all three leptons and the right ones for softest lepton in the event. 
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Figure 8.3: Tracker isolation: a) tt background rejection versus on charged track, 
b) signal acceptance versus p^* on charged track. Events are selected with three leptons of 
Pt > 15 GeV and all of them are required to be isolated. Full and open circles correspond 
to low, L — 10^^ cm~^ s~^, and high, L — 10^^ cm~^ s~^, luminosity regimes, respectively. 
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Figure 8.4: Jet multiplicity and transverse energy distributions for the signal and inter- 
nal SUSY background at mSUGRA point SI (100 GeV, 100 GeV) and the tt background. 
Events are selected with three isolated leptons of px > 15 GeV. 
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Figure 8.5: Jet veto rejection/acceptance versus \rf^^\ in xfx2 searches. Full circles 
correspond to E^:^ > 25 GeV and the open ones to Ej^* > 30 GeV. 
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Figure 8.6: Missing transverse energy distributions in signal and background events 
with three isolated leptons of px > 15 GeV and no jets. 
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Figure 8.7: Dilcpton invariant mass distributions in signal and WZ background events 
with three isolated Icptons of pr > 15 GeV and no jets. 
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9 Possibility of observing h ^ 66 in squark and gluino 
decays 

The usual way to search for the hghtest CP-even Higgs (h) of the MSSM is through its 
direct production by gluon fusion or in association with W, bb or tt and its decay into 77, 
ZZ* — »• 4/, nn or bb JT^, ^ All these production mechanisms are Electroweak. Each 
of these channels has some limitations, such as high tan/5 for the /i/i mode, or high mass, 
rrih ^ 130 GeV and large stop mixing for the ZZ* — > 4/ mode, for example. The most 
general way is to search for inclusive h — 77, which, with Lint = 10^ pb~^, would allow 



us to explore a domain approximately given by > 250 GeV and tan/? > 3 [|I^. With 
the h — i> 77 mode we expect a signal on top of a large irreducible 77 background with a 
signal to background ratio of < 1/20. The instrumental requirements on the ECAL are 
very demanding, a 77 effective mass resolution better than ~1 GeV at rrih — 100 GeV 
p3| . In the associated production modes Wh and tih with h —>■ 77, the S/B ratio is of 
order 1 and calorimeter performance is less demanding: but these channels require > 10^ 



pb [^]. Evaluations of ways to exploit the h — > 66 decay mode in the Wh or tth final 
states leave little hope — no way to get a really significant signal has been found up to 



now 



It is well known that the MSSM h can be abundantly produced in the decays of 
charginos and neutralinos (primarily Xg)- turn, the X2 is a typical decay product of 
squarks and gluinos (see Figs. 2.5a,b) which are also produced abundantly. The idea is 
thus to use strongly produced q/g and to look for the dominant decay mode h — > bb, if 

~^ Xjh is kinematically allowed, and require EJp^'^'^ and jet multiplicity to suppress the 
background |T^, ^ . 

In this study we fix only one of the mSUGRA parameters, = 0, and vary the 
other four: mo, 1711/2, tan/?, sign{fi). The most pronounced variation of physics quantities 
(masses, branching ratios, cross-section) is with mo and mi/2, thus the investigation is 
done in the mo, mi/2 parameter plane for several representative values of tan/3 and sign{fi). 
When M^o — M^o > rrih the main source of Higgs is the decay X2 ~^ Xih- This decay 
has a significant (~20-90%) branching ratio in a large portion of the parameter space, see 
Fig. 2.17. In regions of mi/2 ^ 200-300 GeV (depending on tan/?) the three body decays 
X2 ~^ Xi QQy X2 ~^ Xi ^^^^ take over because the X2 ~^ Xi^ decays are not kinematically 
allowed. The boundary at mo ~ 200-300 GeV, on the other hand, corresponds to the 
opening of two-body cascade modes, e.g. X2 ~^ ^l,r^^ ~^ Xi ■ Figure 2.17 also shows 
that the tan/? = 2, /i < case looks the most promising one, whereas the tan/5 = 10, 30, /i 
< seem the least favorable among the possibilities considered, the X2 ~^ Xi^ branching 
ratio being the smallest on average in these latter cases. Depending on tan/5, there is 
potentially a gap between the upper limit explored up to now or foreseeable at LEP2 and 
Tevatron at mi/2 ~ 100 GeV {ruh ^ 100 GeV) and the lower limit at mi/2 ~ 200-250 
GeV where X2 ~^ Xi^ opens up. This comes about because within mSUGRA there are 
simple mass relations: M^o ~ 2M^o ~ mi/2, and the LSP mass is typically > 100 GeV 
in the regions where X2 ~^ Xi^ is allowed. The range of the h mass varies from ~ 80 to 

120 GeV and the exact upper limit is determined by the order to which the radiative 
corrections have been calculated. 

With increasing mo (at fixed mi/2) squarks become more massive than gluino at some 
point (see, e.g. Fig. 2.16) and thus gluino decays into squarks are not allowed anymore. 
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decreasing significantly the yield of X2 from squark decays. Furthermore, the B{q —)■ 
also decreases with increasing tjiq because the decay q gq plays an increasingly 
important role. These factors, with the variation of the total production cross-section 
(shown in Fig. 2.15) and the variation of the kinematics of the decay chains, require that 
for a quantitative evaluation of the regions where h bb could be observed one needs to 
investigate the parameter space point by point and optimise cuts accordingly |T2| . 

We considered the following SM background processes: QCD 2^2 (including bb), 
tt and Wtb. The pT range of all these processes is subdivided into several intervals 
to facilitate accumulation of statistics in the high-pT- region: 100-200 GeV (except for 
QCD), 200-400 GeV, 400-800 GeV and > 800 GeV. The accumulated statistics for all 
channels correspond to 100 fb~^, except for the QCD jet (instrumental mismeasurement) 
background where the statistics correspond to 0.2 - 50 fb~^ depending on the pT interval. 
It is very time consuming to produce a representative sample of QCD jets in the low-p^ 
range, since the cross section is huge and we need extreme kinematical fluctuations and 
non-Gaussian tail instrumental contributions for this type of background to be within 
the signal selection cuts. On the other hand, our simulation of instrumentally induced 
missing transverse energy is not yet fully reliable as it strongly depends on the still evolving 
estimates of dead areas, volumes due to services, non-Gaussian calorimeter response, etc. 
So we cannot go confidently below ~100 GeV with the cut on E'™'^'', where the QCD jet 
background becomes the dominant contribution. 

Initial requirements for all the SM and SUSY samples are the following: 

• at least 4 jets with E^'* > 20 GeV in | r/^'"* | < 4.5, 

• E^'"' > 100 GeV, 

• Circularity > 0.1. 

No specific requirements are put on leptons. If there are isolated muons or electrons with 
pip^ > 10 GeV within the acceptance they are also recorded. The term "isolated lepton" 
here means satisfying simultaneously the following two requirements: 

• no charged particle with pt > 2 GeV in a cone R = 0.3 around the direction of the 
lepton, 

• T,E^^° in the annular region 0.1 < R < 0.3 around the lepton impact point has to 
be less than 10% of the lepton transverse energy. 

At a later stage we investigate whether selection or vetoing on isolated leptons can be 
used to help suppress the SM or internal SUSY backgrounds. 

In this study we used an evaluation of the 6-tagging performance of CMS from impact 
parameter measurement in the tracker [Q. An example is shown in Fig. 4.11. The 
6-tagging efficiency and mistagging probability for jets is evaluated as a function of Et 
and jet pseudorapidity. It is obtained from a parameterized Gaussian impact parameter 
resolution, dependent on the pixel point resolution, radial position of the pixel layers and 
the effects of multiple scattering on intervening materials, with the addition of a non- 
Gaussian impact parameter measurement tail based on CDF data. Although there are 
some specific assumptions made concerning the tails of the impact parameter distributions 
extending beyond the parameterized Gaussian parts in the presently on-going study 
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of the expected impact parameter resolution, with full pattern recognition and track 
finding in CMS is in a very good agreement with the results used as an input in the 
present study. 

The 6-jet tagging efficiency is parameterized as a function of in 3 intervals of rj: 
- 1, 1 - 1.75 and 1.75 - 2.4. A typical 6-tagging efficiency at ?7 ~ for E^^ = 40 GeV 
is ~ 30%, reaching a maximal efficiency of ~ 60% for high-E'-r jets. Charm-jets have a 
typical tagging efficiency of about 10%, and for light quarks and gluons the 6-tagging 
(mistagging) efficiency reaches a maximum of about 3% for high-ii^T^ jets. At the analysis 
stage the jets produced by the event generator are "6-tagged" according to the tagging 
efficiency for 6-jets and mistagging probability for quark and gluon jets. If not specially 
mentioned otherwise, h — >■ 66 decay jets are tagged only in the barrel \ 1] \ < 1.75 interval 
and in all di-jet mass distributions we take only one pair of 6-jets per event, those closest 
to each other in — (y9 space. No attention is paid to jet "charge" or "leading charge" . 



More details can be found in \12 



As an example let us consider a representative mSUGRA point mo = mi/2 = 500 
GeV, tan/? = 2, yU < 0. Here the masses of relevant sparticles are: rrig = 1224 GeV, m^^ 
= 1170 GeV, m^o = 427 GeV, m^o = 217 GeV, = 89.7 GeV. In Fig. 9.1 we compare 
kinematical distributions for this point with the SM background for Li^t = 100 fb~^. The 
j^mtss (distributions show the most pronounced difference between signal and background, 
this then being the most important variable for the background suppression. A significant 
difference exists also in the jet multiplicity distributions. A specific optimisation of cuts 
on E"™*^*, E^r^^ , etc., is however performed in various regions of mSUGRA parameter space. 

Figure 9.2 shows the h-h (jet-jet) mass distribution for the sum of the SUSY and SM 
events assuming an ideal 6-tagging performance, i.e. 100% 6-tagging efficiency with 0% 
mistagging probability for jets. The cuts are 

SM background shown in Fig. 9.2 separately as a shaded histogram is small compared to 
the internal SUSY background. The bb mass distribution is fitted with a sum of a Gaussian 
and a quadratic polynomial. The width of the signal peak is determined entirely by the 
jet-jet effective mass resolution of the detector (the intrinsic h width varies from 3.2 to 
4.3 MeV over the entire mSUGRA parameter space). The position of the peak is shifted 
to lower mass since energy losses in jets (finite cone size, neutrinos generated in hadronic 
showers, etc.) are not corrected for. The three distributions in the lower half in Fig. 9.2 
are for events in the Mbb mass window 70 - 100 GeV to illustrate the kinematics of the 
h ^ bb decay. The average 6-jet transverse energy is ~110 GeV, and the 6-jets are very 
central as they result from the decay of massive (~1200 GeV) squarks and gluinos which 
are centrally produced. Clearly, 6-tagging beyond | 77-'^* | ~ 1.5 is not very helpful in this 
search. 

If one now applies the nominal 6-tagging performance expectations described previ- 
ously, we obtain the bb mass distribution given in Fig. 9.3a. Under these nominal condi- 
tions and with 100 fb"^ the expected signal significance calculated as S/\fB within a ±lcr 
interval around the peak value is 18.3. The background is dominated by the internal SUSY 
background which is largely irreducible, originating mainly from additional real 6-jets in 
the event, rather than from mistagged jets. An improvement could possibly be obtained 
by taking into account "jet charge" and thus enriching the sample in hb pairs as compared 
to hb + bb/bb. This however requires track reconstruction and charge determination for 
fast tracks in jets, a task still under study in CMS . 

Figures 9.3-9.5 illustrate the expected sensitivity of the h. ^ bb signal to various 
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instrumental factors [|12|. In Fig. 9.3b all 66-combinations per event are included in 
the distribution (instead of just the one with the closest distance in 77 — space) and 
the signal significance is now worse by ~12%. Figures 9.3c and 9.3d show the effect 
of increasing, or decreasing, the 6-tagging efficiency by 15% in absolute value, whilst 
keeping mistagging at the nominal level. The dependence on the 6-tagging efficiency 
is very pronounced, the signal significance varying by respectively ±25%. This is not 
surprising as the 6-tagging efficiency enters quadratically in the number of signal events. 
Figures 9.4a and 9.4b illustrate the effect of the rj acceptance of the silicon pixel detector. 
As can be expected on the basis of Fig. 9.2, increased acceptance beyond | rj^'^^ \ ~ 1.5 does 
not bring a significant improvement as the signal 6-jets are very central. With tagging 
acceptance increasing from | t]^'^^ |= 1.0 to 1.75 and 2.4, signal significance improves from 
17.2 to 18.3 and 18.8, respectively. Figures 9.4c and 9.4d illustrate the dependence on 
tagging purity, where one keeps the nominal 6-tagging efficiency unchanged, but varies 
the mistagging probability for all non-6-jets. Increasing the mistagging probability by a 
factor of 3 degrades signal significance by ~25%. 

Since the observed width of the h — 66 peak is entirely determined by the calorimetric 
jet-jet effective mass resolution, one can expect dependence of the signal observability on 
the energy resolution. This is illustrated in Fig. 9.5. The upper plot is for the nominal 
HCAL performance, corresponding to aE/E = 82%/V^ © 6.5% at = and the lower 
plot corresponds to the HCAL energy resolution degraded to 120%/V^ © 10%. The 
resolution of the Gaussian fit to the Higgs peak is 7.6 GeV in the first case and 11 GeV 
in the second. The difference in S/B between Figs. 9.5a and 9.5b is significant, but not 
dramatic. This is due to the fact that the h. ^ bb pair has a significant boost, ~ 200 
GeV (Fig. 9.2) and the jet-jet opening angle, whose measurement precision is determined 
by calorimeter granularity plays an important role in the effective mass resolution. The 
mass resolution also contains irreducible contributions from final state gluon radiation, 
jet cone size and fragmentation effects. 

Figure 9.6 illustrates the h — > 66 signal observability as a function of uih over the 
allowed range of ~80 to 120 GeV. For this we select the points all at the same mo = wii/2 
= 500 GeV, but different values of tan/3 and /i, with nominal 6-tagging performance and 
Lint = 100 fb^^. One can clearly see that the signal becomes broader with increasing 
Higgs mass since the jet Et > 40 GeV cut-off does not play such a role as for a low-mass. 

It is worth mentioning, that a significant part of the mSUGRA parameter space can 
be explored already at a very low luminosities. For example. Fig. 9.7 shows the expected 
Higgs signal with Li^t = 3 fb^^ for a few points. All the distributions in this Figure are 
obtained with the cuts: EJp^^'^ > 200 GeV and E^^^ > 40 GeV. For comparison, searches 
in the h — 77 channel at the corresponding mSUGRA points would require ~10 times 
higher integrated luminosity [jl2| . 



Let us turn to the main problem, how general is the possibility to observe this h — > 
66 signal in q, g decays. Figure 9.8 shows the domain of parameter space where the h 
bb signal is visible with S/ \^ > 5 for the case tan/5 = 2, = 0, yU < with nominal 
6-tagging performance of CMS, for 10 and 100 fb~^ integrated luminosities. The isomass 
curves for the CP-even (h) and CP-odd (A) Higgses are also shown in Fig. 9.8 by the 
dash-dotted lines. The bold broken line denotes the region where B(x2 Xih) = 50%. 
The shaded regions along the axes denote the present theoretically (TH) or experimentally 
(EX) excluded regions of parameter space not yet including LEP 96/97 results. The LEP2 
and Tevatron (with 1 fb~^) sparticle reaches are also shown by solid lines. With Lint < 1 
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fb"''^ the threshold of visibihty of h ^ 66 at lowest mi/2 = 170-180 GeV corresponds to 
^ 400-450 GeV, i.e. it begins just where the Tevatron searches will stop with ~5 fb^^. 
With 100 fb~^ the reach extends up to rrih — 90 GeV, tua ^ 1500 GeV. Figures 9.9 and 
9.10 show the domains of h signal visibility for tan/3 = 10, < and tan/3 — 30, /i > 0, 
respectively. In both cases there seems to be a significant observability gap between the 
upper reach of LEP2 or the Tevatron and the LHC/CMS low rrih reach using this method 
to look for h. This has to be studied in more detail in the future. 

We conclude that the search for the h — > 66 decay, when the lightest Higgs is produced 
in the cascade decays of the strongly interacting sparticles, seems to be a promising 
channel. The large rejection factor needed to suppress backgrounds and achieve a, S/B oi 
~1 is here provided by the E"^*^* cut. Nothing similar can be obtained in the search for 
the SM Higgs in H ^ 66 nor for MSSM h — > 66 in inclusive h, Wh and tih final states, 
mass in the ~450 to 700 GeV range for this search to be possible. The study carried out 
here in the mSUGRA framework shows that there is a significant domain of parameter 
space, just beyond the g/qmass reach of the Tevatron (~400-450 GeV), where observation 
of the h — > 66 decay would be possible with an integrated luminosity of only 1-3 fb~^. 
This parameter space domain can be significantly extended with 100 fb^^, where g, q with 
masses in 1.5 TeV range are probed. The h mass range from ~80 GeV up to ~125 GeV 
can be covered. 

Our investigations show that for the observation of the h — > 66 signal a calorimetric 
energy resolution of ~ 100%/V^ © 10% is adequate, but it should not be significantly 
worse than this. The signal visibility depends most critically on the 6-tagging efficiency 
and to a lesser extent on the mistagging probability; the acceptance of the 6-tagging 
pixel devices is not a real issue, provided the coverage is no smaller than | r] \ ~ 1.5. 
A significant improvement in 6-tagging efficiency could be obtained with a third pixel 
layer, or by having pixel layer at a radius of 4 cm even for high luminosity running, to be 
replaced every (few) year(s). The mistagging probability depends on the non-Gaussian 
part of the impact parameter measurement distribution, and thus on the overall pattern 
recognition performance in the region close to the beam, e.g. on the balance between Si 
and MSGC layers in the tracker. Optimisation of this aspect might then involve increasing 
the number of Si layers from 4 to 5, possibly 6, reducing correspondingly the MSGC layers 
and deserves a dedicated study. 

If this method is to become a viable alternative to the h ^ 77 search, it is important 
to evaluate how general are the results of the present study, i.e. what happens outside 
the mSUGRA scheme where masses are not so constrained, for example in the MSSM. 
Presumably what is found here in the framework of mSUGRA-MSSM is valid as soon as 
the g/q ^ X2 Xih, or even more generally, the g/q Xi ^ X^h chains are kinematically 
allowed and the XiX^^ couplings non-vanishing. A particular point of interest in this 
respect is the gap at the lower q/g masses, between the domains where SUSY can be 
explored at the Tevatron with 5-6 fb~^, the LEP2 reach in terms of rrih and the lower 
mass reach of this channel, if it exists in mSUGRA, can it be overcome in the MSSM 
where mass relations are less rigid? 
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Figure 9.1: Comparison of mSUGRA signal at mo = mi/2 =500 GeV, tan/3=2, jj, < 0. 
and SM background distributions. 
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Figure 9.2: 6-jet related distributions for the same point as in Fig. 9.1 in case of ideal 
tagging performance. The three distributions in the lower part are for events in the Mf,h 
mass window 70 - 100 GeV. 
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Figure 9.3: Influence of tlie various instrumental factors on signal observability in 
the same parameter space point as in Fig. 9.2 with 100 fb~^: a) with "nominal" b- 
tagging performance and data selection, b) all 66-combinations per event are included in 
histogram, c) 5-tagging efficiency increased by 15% and d) 6-tagging efficiency decreased 
by 15%. The nominal mistagging probability is assumed in all cases. 
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Figure 9.4: Influence of various instrumental factors on h — 66 signal visibility in the 
same parameter space point as in Fig. 9.2: a) tagging acceptance increased compared to 
"nominal" from | r/ | < 1.75 up to | | < 2.4, b) tagging acceptance decreased to | | 
< 1.0, c) mistagging probability increased by a factor of 3 and d) mistagging probability 
decreased by a factor of 2 relative to nominal expectations. 
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Figure 9.5: Effect of the assumed single hadron energy resolution on the signal visibil- 
ity: a) the nominal HCAL performance, as/E = 82%/\/^©6.5% at r/ ~ 0; b) the HCAL 
energy resolution deteriorated to 120%/V^ © 10%. mSUGRA point as in Figs. 9.2-9.4. 
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Figure 9.6: Dependence of the h. ^ bb signal visibility on mh over the allowed range 
~ 80 to 120 GeV, varying tan/5 and sign{fi) at a fixed mo = mi/2 = 500 GeV. Nominal 
CMS instrumental performance, 100 fb~^ integrated luminosity. 



118 



m„ = 200 GeV, = 300 GeV in„ = 300 GeV, m^|^ = 250 GeV 

m (g) = 770 GeV, m (u^) = 705 GeV m (g) = 659 GeV, m (u^) = 642 GeV 




M„ (GeV) 82-4 M., (GeV) 



in„ = 400 GeV, m^,^ = 180 GeV 
m (g) = 505 GeV, m (u,^) = 577 GeV 




A„ = 0, tanp = 2, n<0 



60 I 120 180 



M,JGeV) 



Figure 9.7: Some points in parameter space accessible already with 3 fb ^. Nominal 
6-tagging performance of CMS is assumed. 
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Figure 9.8: 5o" visibility contours of h — > 66 for tan/5 = 2, Aq = and fi < with 10 
and 100 fb^^. See also comments in text. 
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Figure 9.9: 5a visibility contours of h ^ 66 for tan/3 = 10, Aq = and jj < with 10 
and 100 fb~^. See also comments in text. 
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Figure 9.10: 5a visibility contours of h ^ 66 for tan/3 = 30, Aq = and fi < with 
10 and 100 fb~^. See also comments in text. 
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10 Sensitivity to sparticle masses and 
model parameters 

10.1 Establishing the SUSY mass scale 

Once a deviation from the SM expectations is established in a number of final states indi- 
cating the presence of SUSY, the determination of the "SUSY mass scale" would become 
the first objective. The production and decay of high mass strongly interacting particles, 
as would be the case for gluinos and squarks, would lead to hard event kinematics, with a 
number of hard central jets in the final state accompanied by significant missing transverse 
energy due to escaping LSPs and neutrinos. As suggested in we can characterize the 
"hardness" of an event by the scalar sum of transverse energies of the four hardest jets 
and the missing transverse energy: 

= E^ + E^ + E^ + E^ + Ef (10) 

This variable is rather insensitive to the transverse energy flow induced by event pile-up 
since it takes into account only the calorimetric cells which are attached to a jet. 

The following SM processes have been considered as a background to SUSY strong 
production: tt, W/Z + jets and QCD jet production (PYTHIA 2 2 processes including 
bb). Events are selected by requiring: 



• E^''' > max(100 GeV, 0.2E|," 

• > 4 jets with Et > 50 GeV, the hardest with Et > 100 GeV in Ir^l < 5; 

• Lepton veto: events are rejected if they contain an isolated lepton with pr > 10 
GeV. 

Figure 10.1 shows the E^"^ distributions for a representative mSUGRA point (with mo = 
300 GeV, mi/2 = 150 GeV, Aq = -600 GeV, tan/5 = 2, /i > 0) and the SM background. 
The background, which is dominated by tt and W/Z + jets productions, is small compared 
to the SUSY signal. 

The peak value of the E^"^ spectrum for the inclusive SUSY signal provides a good 
estimate of the SUSY mass scale, defined as in lij 



MsusY = mm{Mg, Mg) (11) 

and where Mq is the mass of squarks from the first two generations. Figure 10.2 shows the 
relationship between the peak value of the E^'^ spectrum and Msusy for 100 mSUGRA 
models chosen at random with 100 < mo < 500 GeV, 100 < mi/2 < 500 GeV, -1000 < 
Aq < 1000 GeV, 1 < tan/5 < 12 and sign{^) = ±1 At each point the peak value of 



E^"^ was found by a Gaussian fit. The very strong correlation between the peak value of 
EjV'"^ and Msusy as can be seen in Fig. 10.2 could provide a first estimate of the relevant 
SUSY mass scale. 

10.2 Constraints from ^2 leptonic decays 

Establishing the sparticle spectrum, the measurement of sparticle masses and the model 
parameters will be the next step in uncovering SUSY, once its existence and its mass scale 
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have been established. In i?-parity conserving models undetectable LSP's make this task 
nontrivial. 

As mentioned previously, some specific SUSY decays are particularly useful for the 
determination of sparticle masses. The X2 leptonic decays may not only provide the 
first indication for SUSY in significant regions of mSUGRA parameter space, but also 
allow us to obtain information on some sparticle masses by measuring the endpoint of the 
dilepton invariant mass spectrum and exploiting its specific kinematics. Before discussing 
how to exploit this feature in inclusive two- and three-lepton final states, which is a 
recent development, let us first review the invariant mass spectrum in the exclusive 
XiX2~^ 3/^ + no jets+ EJp^^^ final state, first described in |^ for three-body decays 
X2 Xi- In this case the upper limit of dilepton invariant mass spectrum is: 

MJUT- = M^o - M^o (12) 

Figure 10.3 shows the expected mass spectrum for mSUGRA "Point 1" (see Table 1.1) 
in the exclusive xfx2^ 3/^ +no jets+ (£^™^*) channel arising from direct (Electroweak) 
xfx2 production. The SM and internal SUSY backgrounds are also shown, all for L 
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pb~^ |T^. The edge is at 52 GeV and with these statistics it could be measured with a 
precision better than 1 GeV. This measurement of MJ+f- = M^o — M^o = 52 GeV should 
allow, within mSUGRA, estimation of M^o, M^o, M-±, Mg and mi/2 through eqs. (2) - 
(4) as discussed in section 2.1. 

Even more information can be extracted from the X2 two-body (cascade) decays X2 ~^ 
^^^L_R ~^ ^~^l~Xi- The kinematical upper limit in this case is: 

^(M|o-Mf)(Mf-M|o) 
= ^ — (13) 

i.e. it is sensitive also to the intermediate slepton mass. In this case, however, the 
measurement of only the edge position does not provide information on the masses of all 
the involved particles unambiguously. Nonetheless, as discussed in the following, through 
analyses of some kinematical distributions, the masses of Xi? X2 of the slepton can 
be determined with reasonable precision [Q. Therefore, with two-body decays one can 
constrain mi/2 and niQ through eqs. (2), (3) and (5)/(6), if the information on tan/? is 
available, for example from the Higgs sector of the model. 

Direct leptonic decays of the X2 dominate below mi/2 ^ 200 GeV for nearly all values 
of mo; there is also a small region of high mi/2 where direct three-body decays are open, 
as can be seen from Fig. 2.3a. Cascade leptonic decays of X2 (Pig. 2.3b) or 

(Fig. 2.3c) occur at almost any mi/2 as soon as mo ^ 0.5mi/2, where X2 is heavier than 
the sleptons. 

The expected position of the edge as a function of mi/2 is given in Fig. 10.4 for two 
different values of the common scalar mass: mo = 400 GeV and 100 GeV. The box size in 
the Figure is proportional to the branching ratio B{x2 + invisible). At mo = 400 

GeV the X2 only direct three-body decays for all values of mi/2 up to mi/2 ^ 180 GeV, 
where "spoiler" modes open up. Measurement of thus yields the common gaugino 

mass parameter mi/2 unambiguously in this region. For mo = 100 GeV the situation 
is more complicated. Here at different values of mi/2 different leptonic decay modes of 
X2 dominate: direct decays up to mi/2 ^ 130 GeV; cascade decays via 1^ at 130 GeV 
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~ "^1/2 ~ 250 GeV and, finally, cascade decays via li at mi/2 ^ 250 GeV. In some regions 
two X2 decay channels coexist. 

The essential issue, for the exploitation of an observed edge, is knowing whether the X2 
decay chain is three-body or two-body. In general, the pt spectra of leptons from X2 two- 
body decays are more asymmetric compared to the three-body decays. To characterize 
the asymmetry we introduce the variable 0: 

A = ^ ^ (14) 

Pt "t Pt 

where p^"-^ (p™") corresponds to the lepton of maximum (minimum) transverse momen- 
tum. 

An example of decay type determination is discussed for the mSUGRA point rriQ = 100 
GeV, mi/2 = 150 GeV. The sparticle masses are M^o = 135 GeV, M^o = 65 GeV M^^ = 

120 GeV and X2 decays via X2 ^ ~^ ^^^^ + Xi with a probability of 0.54. Figure 10.5 
shows the dilepton mass spectrum for this parameter space point in the inclusive 3 lepton 
events for an integrated luminosity of Lint = 5 x 10^ pb~^. The expected SM background 
is also shown. Of the three leptons with px > 15 GeV only the two of opposite charge and 
same flavor, which enter the invariant mass distribution, are required to be isolated. A 
very spectacular edge is situated at ~ 52 GeV. There is a second and much weaker edge 
at ~ 69 GeV which is due to direct three-body decays with B(x2 ~^ l^l^ + Xi) = 0.05. 
Thus at this mSUGRA point there are two leptonic decay modes with very different 
branching ratios. To identify the decay chain responsible for the first edge, we look at 
the asymmetry A distribution of the lepton pairs with Mi+i- < 52 GeV. Figure 10.6 
(full line) shows the A distribution for these events. The pronounced asymmetry in pt 
indicates the cascade nature of decays. If we now pick up the lepton pairs from the 
mass interval 55 < Mi+i- < 69 GeV below the second edge, then the corresponding A 
distribution peaks near zero, as seen in Fig. 10.6 (dashed-dotted line), indicating the 
three-body decay type of these events. To illustrate the generality of this approach, two 
more examples of A distributions are also shown in Fig. 10.6 for mSUGRA parameter 
space points mo = 100 GeV, mi/2 = 400 GeV with pure 2-body decays (dashed line) and 
mo = 400 GeV, mi/2 = 150 GeV with pure 3-body decays (dotted line), respectively. 

Once the decay type is determined, the next step is to extract the masses of the 
particles involved, the X25 Xi cind /. In our example case of two-body decay the l^l^ 
kinematical upper limit is 

^(M|o-Mf)(Mf-M|o) 
MTff = ^ ^ ^ = 52 GeV. (15) 

This equation can be satisfied by an infinite number of M^o, M^o and mass combi- 
nations. To find a solution, we assume M^o = 2M^o and generate xH ~^ ^^^^ ^ ^^^~Xi 
two-body sequential decays for various assumptions on Mj (M^o), e.g. from 70 GeV to 150 
GeV in 10 GeV steps, with the Xi mass constrained to provide the "observed" position 
of the edge. For each value of Mj, one obtains two series of solutions for equation (15). 
Some of them are shown in Fig 10.7 in terms of dilepton invariant mass and p^^-asymmetry 
distributions, in Figs. 10.7a,c for "low mass" solutions and Figs. 10.7b,d for "high mass" 
solutions M. The solid lines in Figs. 10.7a-c correspond to the "observed" spectra and the 
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dotted histograms are the results obtained from eq. (15). Clearly, the "observed" 
invariant mass spectrum itself and the p^-asymmetry distributions allow us to eliminate 
the "high mass" solutions of eq. (15). To find the best combination of M^o, M^o and 
Mj among the "low mass" solutions, we perform a x^-test on the p^-asymmetry distri- 
butions, taking into account only the difference in shapes, but not the normalization. 
The result is shown in Fig. 10.8. The horizontal line in this Figure corresponds to the 
expected uncertainty of simulations due to the detector resolution, background estimates, 
initial/final state radiations, etc. We obtain the following precisions on masses: SM^o < 
5 GeV, 6Mi < 10 GeV. The use of the total number of observed events, as well as the 
dilepton invariant mass spectrum itself in a combined x^-test would further improve these 
results. 

As shown above for the mSUGRA point under consideration, the second edge in 
the dilepton mass spectrum at MJ+p = 69 GeV (see Fig. 10.6) is due to direct three- 
body decays of X2- Using now two measured edge position values and assuming again 
M^o = 2M^o, two solutions for the slepton mass can be directly obtained Mj = 120 GeV 
and Mj = 77 GeV. The corresponding dilepton mass spectra for these two solutions are 
shown in Fig. 10.9. The clear difference in the predicted spectra (the first edge in Fig. 10.5 
which proceeds through / vs. Fig. 10.9) allows us to eliminate the Mj = 77 GeV solution. 
At this particular mSUGRA point the use of two observable edges provides a precision of 
SM^oj< 1 GeV il|,i. 

Another example of a double edge is given in Fig. 10.10 for the mSUGRA point 
(50 GeV, 125 GeV). Here the sparticle masses are M^o = 116 GeV, M^o = 55 GeV, 
= 110 GeV, = 78 GeV. In this case two two-body decays, via left and right 
sleptons coexist with the comparable branching ratios B{x2 ^^^J ^~^^~Xi) = 0.037 
and B(x2 ^ ^^^r ^^^~Xi) = 0.013 and their analysis could proceed as indicated above 
providing a strong constraint on the underlying model. 

10.3 Determination of the squark mass 

The observation of an edge in the dilepton mass spectrum, resulting from X2 leptonic 
decays allows not only the determination of x?) X2' (0 masses, but also enables the 
momenta to be fixed. At the upper limit of the l^l^ spectrum in direct three-body decays 
X2 ~^ ^~^^~Xij the Xi is produced at rest in the X2 ^^st frame. Assuming knowledge of the 
Xi mass (or, equivalently, knowledge of a relation between Xi and X2 masses), one can 
then, for events at the edge, reconstruct the X2 momentum vector in the laboratory 
frame: 

p^o = il + M^o/Mi+i-)xpi+i- (16) 

Once the four-momentum vector of the X2 is determined, one can search for a resonance 
structure in, e.g. the X2 + j6t(s) invariant mass distributions. This technique can be 
put to a good use to reconstruct g/q masses [|8| as the next-to-lightest neutralinos are 
abundantly produced in gluino and squark decays, see e.g. Figs. 2.5a,b. 

An example of the q mass reconstruction is given for mSUGRA point with itlq = 300 
GeV, mi/2 = 150 GeV, Aq = —600 GeV, tan/? = 2 and /i > 0, discussed in section 10.1. 
The masses of relevant sparticles are 463 GeV, 495 GeV, 151 GeV, 425 GeV, 112 GeV and 
57 GeV for qi^R (squarks of the first two generations), ^2, ^1, 9, X2 Xi^ respectively, and 
the total SUSY cross-section (452 pb) is largely dominated by c^t^i^ — 195 pb, f^g^^ = 109 
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pb, agg = 77 pb productions. Gluinos are lighter than squarks of the first two generations, 
but their decays to stop and sbottom are kinematically allowed, giving several 6-jets at 
the end of the decay chain. Right squarks qr predominantly decay to Xi + Q or to g+q. 
Finally, the Left squarks qi have the interesting for us decay mode to X2+ 1 with a 
branching ratio of 0.27. To extract this latter decay chain from internal SUSY and SM 



backgrounds we have adopted the following procedure p9 



• Require two isolated leptons withpr > 10 GeV; and an angular separation ARi+i- = 
VAr/^ + A(p'^ < 1. 

• The scalar sum of calorimeter cell transverse energies and of muon transverse mo- 
menta, E£^|?'° +Pt > 600 GeV to ensure production of hard events. Require E^'^'^'^ 
> 100 GeV. 

• The invariant mass has been reconstructed and the edge position at ^54 GeV 
has been observed. The p^-asymmetry distribution of dilepton pairs indicates it is 
due to a the three-body decay of type X2 ~^ ^^^~Xi- 

• Events with a dilepton mass close to the kinematical limit, i.e. 49 GeV < Mi+i- < 54 
GeV, have been selected for further analysis; the X2 momentum has been determined 
using eq. (16). 

• The invariant mass of the X2 each non-6-jet with Ej^'^ > 300 GeV has then 
been reconstructed. Such a hard jet is unlikely to originate from a QCD radiation. 
Furthermore, due to the large mass difference Mg^ — M^o = 351 GeV the Et- 
spectrum of quark-jets from the qi decay is harder than both the inclusive SUSY 
and SM jet spectra. Hence this requirement enhances the q^ signal. 

Figure 10.11 shows the reconstructed X2 + J^t mass spectrum for SUSY and SM events 
for an integrated luminosity of L^nt = 10'^ pb~^. As a source of the SM backgrounds we 
have considered ti, W/Z -|- jets, WW, ZZ, ZW and bb processes. Their contribution is 
neghgible compared to the internal SUSY background. A peak at ~450 GeV is clearly 
visible and a Gaussian plus polynomial fit yields a value of 447 ± 4 GeV (statistical error) 
for the reconstructed q^ mass. Significant systematic uncertainties are associated with 
this measurement due to the jet reconstruction algorithm, energy scale, approximations 
in determining the X2 four-momentum, in particular the dilepton mass interval which 
has been chosen for analysis, etc. The overall systematic uncertainty is estimated to be 
about ±5%. A much more detailed analysis is needed to understand fully the achievable 
precision 
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As a concluding remark, the observation of a resonance peak in the X2 + J^^ invariant 
mass spectrum would be interpreted as a squark undergoing a two-body decay to X2 
a quark. The technique described here is also applicable for gluino mass reconstruction, 
for example in ^ — qq'x2 with subsequent decay of X2 ^ Xi- 



10.4 Sensitivity to model parameters in 3/ + no jets + Ej^'^^^ final 
states 

Since information on mi/2 can be obtained from the edge position in regions of 
parameter space where it is visible, one can then attempt to constrain the mo parameter 
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Table 10.1: SUSY cross-sections (in fh) after selection. Events are counted below 

"observed" edge at 54 GeV. 



Point 


Xl X2 


99 9Q QQ 


9X QX 


XX 


// lu vv 




( i.O 


^ A ^ 


i i .o 


K Q 

o.y 


1 /I 


(150,100) 


53.8 


15.9 


5.3 


2.6 


5.3 


(200,100) 


45.6 


11.9 


6.0 


4.0 


0.0 


(300,100) 


29.7 


0.5 


1.0 


0.5 


0.0 


(400,100) 


23.0 


0.0 


0.0 


0.0 


0.0 


(1000,100) 


15.7 


0.0 


0.0 


0.0 


0.0 


(1800,100) 


12.3 


0.0 


0.0 


0.0 


0.0 



Table 10.2: SM cross-sections (in fb) after selection. Events are counted below 

"observed" edge at 54 GeV. 



wz 


ZZ 


tt 


Wt6 




bb 


1.42 


0.48 


1.07 


0.46 


0.27 


< 0.15 



from the event rates. Let us first discuss the case of 3/ + no jets + E^^^^ events |]T0|- For 
a fixed mi/2 the a-B{xtx2 — ^ 3/ + invisible) decreases monotonically with increasing 
in the region where X2 ^as direct three-body decays. Table 10.1 gives the expected SUSY 
event rates for a series of points with different mo but with fixed mi/2 = 100 GeV, i.e. 
corresponding to approximately the same M™^f:'. The remaining parameters are tan/5 = 2, 
Aq = 0, fi < 0. All expected SUSY sources of such events are included. The apphed cuts 
are (Set 1 in section 8): 

• 3/ with p'j, > 15 GeV in \ri''\ < 2.4(2.5) for muons (electrons); 

• Isolation: no track with px > 1.5 GeV in a cone R = 0.3 about the lepton direction, 
for all leptons; 

• Jet veto: no jet with E^""* > 25 GeV in \r]^''^\ < 3.5; 

• Z-mass cut: no lepton pair with Mi+i- > 81 GeV. 



Table 10.3: statistical {Lint = 10^ pb ^) and systematic errors on mo measurement for 

mi/2 ~ 100 GeV. 



Source 


mo G [100 GeV; 300 GeV] 


mo G [300 GeV; 400 GeV] 


statistical error 


8 ^ 20 GeV 


20 ^ 55 GeV 


edge measurement precision 


4^8 GeV 


8 ^ 25 GeV 


uncertainty on SM background 


2.5 ^ 9 GeV 


9 ^ 30 GeV 


uncertainty on luminosity 


27 ^ 35 GeV 


35 ^ 85 GeV 


statistical © systematic error 


30 ^ 45 GeV 


45 ^ 105 GeV 
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Table 10.4: SUSY cross-sections (in fb) after selection. Events are counted below the 

"observed" edge at 72 GeV. 



Point 


Al A2 


00 OO OO 
iJiJ iJH HH 


OY f7Y 


YY 
A A 


// Iv w 


(150,150) 


26.3 


0.3 


1.5 


1.3 


2.9 


(200,150) 


19.9 


0.7 


0.8 


0.8 


0.8 


(300,150) 


10.7 


0.5 


0.1 


0.1 


0.0 


(400,150) 


6.4 


0.0 


0.0 


0.1 


0.0 


(1000,150) 


3.6 


0.0 


0.0 


0.0 


0.0 


(1800,150) 


3.3 


0.0 


0.0 


0.0 


0.0 



Table 10.5: SM cross-sections (in fb) after selection. Events are counted below the 

"observed" edge at 72 GeV. 



wz 


zz 


tt 


Wt6 


Z66 


hh 


1.53 


0.42 


1.66 


0.50 


0.33 


< 0.15 



Prom the selected events we count the ones with Mi+i- < MJ+p — 54 GeV. Con- 
tributions from different SUSY sources - X1X2 direct production, strong production, 
gluino/squark associated production, gaugino pair and slepton pair productions - are 
given separately. The main contribution is always due to xf X2 direct production, but 
below mo ~ 300 GeV other SUSY processes also contribute significantly. 

The expected SM event rates after selection cuts arc given in Table 10.2. All SM pro- 
cesses together amount 3.85 fb, which is always smaller than the X1X2 direct production 
rate at mi/2 = 100 GeV, and much smaller than that for rriQ < 300 GcV (Table 10.1). 

Tabic 10.1 shows that the data selection does not wash out significantly the event 
rate sensitivity to mo expected from the behavior of cr-B for xtx2 direct production 
(Fig. 2.4). Figure 10.12a shows Nsusy versus mo for mi/2 = 100 GeV, where Nsusy 
is the SM-background-subtracted number of 31 + no jets + events surviving cuts 

for an integrated luminosity of Lint = 10^p6~^. The error band includes both statistical 
and systematic uncertainties; The systematic uncertainties include effects due to the edge 
measurement precision, the uncertainty on the rate of SM processes and a 10% uncertainty 



Table 10.6: statistical {Lint — 10^ pb ^) and systematic errors on mo measurement for 

mi/2 ~ 150 GeV. 



Source 


mo e [150 GeV; 300 GeV] 


mo e [300 GeV; 400 GeV] 


statistical error 


5 ^ 10 GeV 


10 - 


- 18 GeV 


edge nieasureniont precision 


2^8 GoV 


5 - 


- 8 GoV 


uncertainty on SM background 


9 24 GeV 


24 - 


- 55 GeV 


uncertainty on luminosity 


28 30 GcV 


30 - 


- 40 GcV 


statistical © systematic error 


32 ^ 40 GeV 


40 - 


- 75 GeV 
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on the luminosity measurement. The resulting error on mo is 30 to 45 GeV for 100 GeV 
< mo < 300 GeV and 45 to 105 GeV for 300 GeV < mo < 400 GeV. Above mg ~ 400 
GeV sensitivity to mo is lost (see also Fig. 2.4). Table 10.3 gives the statistical and the 
various systematic uncertainty contributions to the mo measurement error; the error on 



the integrated luminosity is the dominant contribution |10 . 

For mi/2 = 150 GeV the expected event rate is much lower than for mi/2 = 100 
GeV, and Lint = 10^ pb~^ is needed to obtain a reasonable precision on mo. The cuts 
applied are the same as for the mi/2 = 100 GeV points, except for a relaxed jet veto 
requirement (Set 2 in section 8): no jet with E^* > 30 GeV in \ri^'^^\ < 3. Selected events 
are counted below the edge at 72 GeV. The event rates for the mSUGRA points and 
the SM background processes are given in Tables 10.4 and 10.5. After cuts the sum of SM 
processes is 4.6 fb. The SUSY/SM ratio exceeds 1.4 for mo ^ 400 GeV. Figure 10.12b 
shows NsusY versus mo for mi/2 = 150 GeV assuming Lint = I0^pb~^. The sensitivity to 
mo is 32 to 40 GeV for 150 GeV < mo < 300 GeV and 40 to 75 GeV for 300 GeV < mo < 
400 GeV. Contributions from the different sources of uncertainty are given in Table 10.6 
The error on the luminosity is still the main source of uncertainty, although the error on 
the SM background now has a significant influence due to the lower signal/background 
ratio compared to the mo=100 GeV case. 



10.5 Sensitivity to model parameters in inclusive /+/ + E"™** + 
jets final states 

More generally, information on event rates from the observation of an edge in inclusive 
studies can also be used to constraint the model parameters. Let us discuss the case 
of Z+Z" + E™'*'' (+ jets) final states 0. Within mSUGRA the expected edge position 
]\/[rnax ^^iq dileptou mass distribution can be obtained from equations (12), (13). Figure 
10.13 shows the contours of expected values of M^'^m in the (mo,mi/2) parameter plane. 
Different lines with same values of MJI^p belong to domains I, II and III which correspond 
to the three possible decay modes of X2 l^l~Xi ^^sl states. The regions of MJ+p 
accessible at LHC are: 

for X2 ^ Xil^l' - 50 GeV < Ml^f^ < 90 GeV (I) 

X'i^ld - MJ^I^ZlOGeV (II) 

hi - MJI^p > 20 GeV (III) 



The first case is limited by the appearance of the spoiler modes X2 ~^ h,^{Z'^)Xi- In 
the last two cases the upper limit on accessible M™"j? is determined only by the available 
statistics. A measurement of in the inclusive dilepton mass distribution, with a 

single edge, thus constrains the model parameters in general to three lines in the (mo, mi/2) 
parameter plane. In the case of > 90 GeV the constraint is stronger, as there are 

just two possible lines. The most favorable case is when the measured MJ+'^T: value is 
large, MJ+^H: > 180 GeV. Then one is left with a single line in the (mo,mi/2) parameter 
plane. The observation of two edges due to the simultaneous presence of two types of X2 
decay modes to l~^l~Xi would give even stronger constraints, see Refs. [[?|, P, p| . 



The proper line in the (mo, 1TI1/2) plane can be determined by the method discussed in 
0. When the correct M]+1- line is identified, the next step is to find the point (mo, mi/2) 
on this line. In general the cross section decreases with increasing mi/2 and mo, thus 
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Table 10.7: values (in GeV) at the investigated (mo,?72,i/2) points from domain I, 

xl-^x\ + i^ + i-. 





(120,160) (130,160) (180,160) (200,160) (220,160) (240,160) (290,160) (350,160) 


j^rnax 


74 


74 


74 


74 


74 


74 


74 


73.7 



Table 10.8: MJ^p values (in GeV) at the investigated {rriQ, 1711/2) points from domain II, 



(mo,mi/2) 


(80,162) 


(90,170) 


(105,180) 


(110,187) 


(120,195) 




74 


74 


73 


75 


73 



Table 10.9: MJ+p values (in GeV) at the investigated {mo,mi/2) points from domain 

III, xl^it + i^^x^ + i^ + i~- 



(mo,mi/2) 


(20,195) 


(40,210) 


(60,230) 


(80,255) 




73 


73 


73 


73 



the study of the event rate along the corresponding M™"5 lines will determine the point 
in parameter space. As an example of the general situation let us discuss the case of 
]\^rnax = 74 -|- X GcV, with three lines corresponding to domains I, II and III, respectively. 
The (mo,mi/2) points analysed are given in tables 10.7-10.9. 

As this is an inclusive analysis, we assume a low luminosity of Lint = 10^ pb^^. We 
discuss first domain III, where the situation is simplest. To minimize the uncertainties 
due to background, we require P'^'^ > 15 GeV and E^'^^^ > 130 GeV. The dependence of 
the resulting event rate on mo is shown in Fig. 10.14. The errors are calculated taking into 
account the statistical error and assuming a systematic error of 30 % for the background 
uncertainty; a systematic error due to the precision of the edge position measurement 
is also taken into account. From the observed event rate mo can be determined with 
good precision, 6mo ~ 4 GeV. The parameter mi/2 is then given by the M™"f^-line in the 
(mo,mi/2) plane (Fig. 10.13). The precision obtained in such a way is Smi/2 — 4 GeV. 

In domain II, the event rate along a line of definite M™"5 first increases and then 
decreases with increasing mo. Fig. 10.15. This is mainly due to the change in the branching 
ratios. The event rate thus does not determine mo uniquely: a given event rate corresponds 
in general to two mo values. The ambiguity can, however, be solved at high luminosity 
Lint = 10^ pb~^, when two edges in the Mi+i- distribution can be observed 0. 

For domain I, the mo dependence of the event rate is shown in Fig. 10.16a, again for 
^max ^ 74 ± 1 GeV. There is a steep increase of the rate at mo ^ 120 - 130 GeV. This 
is due to the decay channel X2 ~^ ^^^ Xi just opening up in this region. As can be seen 
from Fig. 10.16a, there is an ambiguity in the determination of mo if the event rate is in 
the region 3700 < Nev ^ 5600 or 120 GeV < mo < 240 GeV. Additional information 
can however be obtained from the average number of jets <Njet> in the events. Figure 
10.16b shows <Njf.t> as a function of mo. <Njet> is increasing with mo as more jets are 
produced with increasing squark-mass. With the measured <Njet> we can resolve the 
ambiguity in the region 120 GeV < mo ^ 240 GeV and determine mo with 6mo ^7-3 GeV. 
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Figure 10.1: E^'^ distribution for both the inclusive SUSY and the SM backgrounds after 
the event selection cuts have been applied. 




Figure 10.2: The relationship between the peak value of the E^'^ distribution and the 
SUSY mass scale, as defined in the text. 
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mass determination in 3-^- + no jets + E[^'ss 
final state from %^ production 
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Figure 10.3: Dilepton invariant masses distribution for mSUGRA point (mo=200 GeV, 
mi/2=100 GeV) with = iC pb"^ in the 3/ + no jets + events. Contributions 

from SM and SUSY backgrounds are also shown. 
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Figure 10.4: Correlation between MJ^^- and mi/2 for fixed mo = 400 GeV and 
GeV. The otfier mSUGRA parameters are: tan/3=2, Aq — and // < 0. 
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mo = 1 00 GeV, mi/2 = 150 GeV, tanp = 2, Aq = 0, ^< 
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Figure 10.5: /"*"/" mass distribution for mSUGRA point (mg = 100 GcV, mi/2 = 150 
GeV); the other mSUGRA parameters are: tan/3==2, = and /j, < 0. The shaded 
histogram corresponds to the SM background. 
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Figure 10.6: Lcpton transverse momentum asymmetry distributions in decays. 
mSUGRA parameters (mo,mi/2) are the following: 100 GeV, 150 GeV (full and dashed- 
dotted lines); 100 GeV, 400 GeV (dashed hue) and 400 GeV, 150 GeV (dotted line). 
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Figure 10.7: The dotted lines correspond to the invariant mass spectra and transverse 
momentum asymmetry distributions in X2 two-body decays for several "low mass" and 
"high mass" solutions of eq. (10). The full line corresponds to the initial (observed) 
spectrum. 
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Figure 10.8: Results of the test of the shape of pT-asymmetry distributions as a 
function of slepton mass. The corresponding values of M^o are also indicated. 
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Figure 10.9: Predicted invariant mass spectra for two values of slepton masses: 
Ml = 120 GeV (full line) and = 77 GeV (dashed line); mSUGRA parameters as in 
Fig. 10.6. 
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Figure 10.10: Expected mass spectrum for mSUGRA Point: rriQ = 50 GeV, 

mi/2 = 125 GeV, tan/5=2, Aq = and < 0. a) Leptons with pt > 15 GeV in |?7| < 2.4 
are considered. Two leptons which enter invariant mass spectrum are required to be 
isolated; b) as in previous case, but all three leptons are isolated. The dotted histogram 
corresponds to the SM background. 
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Figure 10.12: a) Number of SUSY 3/ + no jets + {El^''') events with U^t = 10 
as a function of common scalar mass mo for mi/2 ~ 100 GeV. b) Number of SUSY 
31 + no jets + (E^'^^^) events with Lint — lO^pb"^ as a function of common scalar mass 
mo for mi/2 ~ 150 GeV. 
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Figure 10.13: Contour lines of expected MJ+I- values (in GeV) in the invariant dilepton 
mass distribution corresponding to the three different decay modes in the region of the 
(mo,mi/2) parameter plane accessible with 10^ pb~^. 
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Figure 10.14: Expected 2l + E^''' event rate with = 10^ pb^^ along MJ^p = 74±1 
GeV contour line in domain III (x2 ~^ ^l^^ ~^ Xi^^^ ) as a function of mo- Event selection 
criteria are: pj'' > 15 GeV, EJp''' > 130 GeV and Mi+i- < M^I^p at corresponding points. 
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Figure 10.15: Expected 2/ + £'™** event rate with L 
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Figure 10.16: a) Expected 21 + E^''' event rate with L^nt = 10^ pb"^ along MJ^p = 
74 ± 1 GeV contour line in domain I Xi^~^^~) ^ function of mo. Event selection 

criteria are: prj^'^ > 15 GcV, i?™''* > 130 GeV and M;+;- < MJ+p at corresponding 
points; b) Average number of jets (Ej^^ > 30 GeV, | rjjet | < 3) at investigated points from 
domain I. The numbers in parenthesis show the masses of squarks and gluinos (in GeV) 
at corresponding points. 
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11 Summary and conclusions 



We have investigated the abihty of the CMS detector at the Large Hadron CoUider to 
discover and characterise supersymmetry. As a benchmark model, we used the minimal 
supergravity-inspired supersymmctric standard model (mSUGRA). In the investigated 
scenarios the lightest supersymmetric particle (LSP) is stable and is the Xi- 

Discovery of supersymmetry at the LHC will be relatively straightforward. It would 
be revealed by large excesses of events over the standard model expectations in a number 
of characteristic signatures - for example missing Et plus jets, with one or more isolated 
leptons; an excess of trilepton events; a very characteristic signature in the l^l" invariant 
mass distribution. 

Because of the importance of missing Ex signatures, considerable effort has been 
expended in optimising the performance of the CMS detector for EJp^'"^^ measurements. 
In particular we have tried to minimize sources of non-Gaussian tails in the calorimeter 
response due to cracks and dead materials. 

For the cases we have studied in detail and for 10^ pb~^ integrated luminosity, 

• Squarks and gluinos can be discovered up to masses well in excess of 2 TeV, i.e. 
for a parameter space region mo ^ 2 — 3 TeV and mi/2 ^ 1 TeV (see Fig. 11.1). 
This covers the entire parameter space over which supersymmetry can plausibly be 
relevant to electroweak symmetry breaking without excessive fine-tuning. Over the 
lifetime of LHC experiments these searches would probe the ~ 3 TeV region. 

• Charginos and neutralinos can be discovered from an excess of events in dilepton 
or trilepton final states. Inclusive searches could give early indications from their 
copious production in squark and gluino cascade decays just from the /"*"/" spectrum 
shapes. With 10^ pb~^ the dilepton edge due to the X2 leptonic decays can be seen 
up to mi/2 ~ 900 GeV for mo ^ 400 GeV (see Fig. 11.2). Applying isolation 
requirements and a jet veto would enable direct chargino/ncutralino production to 
be selected; the mass reach in this case covers up to ~ 170 GeV for xf, X2 

~ 90 GeV for Xi- The region of sensitivity (Fig. 11.1) is entirely determined by the 
production cross section and decay branching ratios rather than by any detector 
limitations. 

• Directly-produced sleptons can be detected up to masses of ~ 350 GeV using dilep- 
ton signatures. Again, the use of a jet veto is crucial, and the lepton opening 
angle can be used to discriminate between the slepton signal and other SUSY back- 
grounds. The region of sensitivity (Fig. 11.1) is once again largely determined by 
the production cross section and decay branching ratios. It is also interesting to 
remark that the domain of parameter space explorable through direct slepton pro- 
duction searches matches closely the x^ cosmological dark matter favored region 
of parameter space. This is due to rehc Xi densities being largely determined by 
slepton masses. Observation of a dilepton edge up to mi/2 ~ 900 GeV gives indirect 
access to slepton production for masses up to ~ 700 GeV. 

• Squark and gluino production may in addition be a copious source of Higgs bosons 
through the cascade decays of the SUSY intermediate particles. In this case it will 
be relatively simple to reconstruct the lightest SUSY Higgs from its decay h bb 
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by requiring two 6-tagged jets. Unlike the case of the Standard Model Higgs in this 
decay mode, good signal to background ratios of order 1 can be obtained selecting 
SUSY events using a hard cut on EJp^'^^ thanks to the escaping LSP's (Fig. 11.2). 

Once super symmetry has been discovered, the challenge will be to determine the 
sparticle spectrum and to pin down model parameters. Over much of the parameter space, 
we have shown that, within this model, neutralino masses can be determined using the 
dilepton invariant mass distribution. Details of the X2 decay scheme may be extracted 
from the lepton p^^-asymmetry, for example. What can exactly be measured and with 
what precision depends on the region in mo,mi/2 space Nature has chosen, and what 
topology/final state we use. As visible from Fig. 11.1 some regions of parameter space 
can be explored through several channels allowing for example simultaneous observation 
of g, g and of xf , X2 ^'^d I direct production giving thus access to a larger spectrum and 
providing more constraints on the supersymmetric models. 

One of the attractive features of i?-parity conserving SUSY models is that the x? ^ 
lightest supersymmetric particle is a plausible candidate for the cosmological dark matter. 
The region of parameter space where this is true for tan/3 = 2 is shown in Fig. 11.1; as 
already mentioned, it is covered by both the squark and gluino searches and the slep- 
ton search, and overlaps the region where inclusive charginos and neutralino signals are 
visible (Fig. 11.2). Specifically, in q/g searches even with 10^ pb~^ the x\ (WIMP) is 
probed up to ~ 350 GeV for any mo, up to ~ 450 GeV for mo ^ 400 GeV, i.e. well 
beyond the cosmologically plausible domain. The cosmologically favored region (hatched 
area in Fig. 11.1) is embedded deeply within the explorable domain of parameter space, 
implying that, if SUSY is responsible for cosmological dark matter, this could hardly es- 
cape discovery at the LHC, at least for tan/5 values investigated. This result is probably 
valid much more generally than implied by this specific model since there is a large safety 
margin, and a variety of experimental signatures/final states with different sensitivities, 
backgrounds and systematics. In this respect a detailed investigation of large tan/? values 
would certainly be of interest. 

We have not yet explored other supersymmetry scenarios in detail, but some general 
comments may be made. Generator-level studies 0] have shown that even if the lightest 



neutralino decays hadronically, the missing Et plus jets and isolated leptons signature 
will remain sensitive to squark and gluino production. We also expect that the inclusive 
multilepton signature would be robust against changes in the SUSY decays. Gauge- 
mediated SUSY breaking models, which predict decays of the lightest neutralino into a 
photon plus a gravitino, would give in our detector striking signatures of diphotons plus 
missing Et that would be easily detected . 



We therefore believe that it will be very hard for electroweak-scale supersymmetry, if 
it exists, to escape detection at the LHC. If supersymmetry is still not found after the 
LHC has accumulated few times 10^ pb^^ of data then it is essentially excluded at the 
electroweak scale, and some other new physics is presumably at work. On the other hand 
if supersymmetry is discovered at LEP or the Tevatron, there will be a rich program of 
SUSY studies to be carried out at the LHC since so much more of the Higgs and sparticle 
spectrum will be accessible. 
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Expected reach in various channels 
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Figure 11.1: 5a significance contours at L^nt = 10^ pb ^ for different SUSY cliannels. 
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mSUGRA parameters: tanp= 2, Ao= 0, |li<0 
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Figure 11.2: mSUGRA parameter reach in h — > 6fc decays at 10^ pb~^ and 10^ pb~^ 
integrated luminosities. The reach in via observation of edge(s) in dilepton mass 
distribution is also shown. 
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